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ABSIRUCT 


\ 


There  existed  a  leed  for  an  interactive  program  chat 
would  provide  the  usar  assistanca  ir  solving  applications  of 
linear  control  theory.  The  Liaaar  control  Program  (LIMCON) 
and  its  user's  guida  satisfy  this  need.  A  series  of  ter. 
interactive  programs  are  presented  which  permit  the  user  to 
carry  out  analysis,  design  and  simulation  of  a  broad  class 
of  linear  ccntrol  problems. 

LISCON  consists  of  two  groups:  matrix  manipulation, 
transfer  function  and  time  response  programs;  and  modern 
controls  programs.  Examples  for  each  are  worked  within  each 

terminal  session  section. 
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iNIgDOgCTION 


The  purpose  of  this  thesis  was  to  update  an  existing 
program  which  provides  assistance  in  solving  computational 
problems  associated  with  the  study  and  application  of  linear 
control  theory,  Tha  Linear  Controls  Program  (LINCOm  was 
first  developed  by  Melsa  [  1]  and  adapted  for  batch  use  at 
NPS  by  Desjardins  [2],  Although  the  original  intent  of 
this  thesis  was  simply  to  take  Desjardins'  adapted  version 
of  Melsa* s  LINCON  and  further  adapt  it  by  making  it 
interactive,  LINCON  soon  began  to  grow  as  other  routines 
were  incorporated,  as  will  be  noted,  until  its  present  form 
was  achieved. 

LINCOM,  as  such,  is  a  higi  level  applications  software 
system  made  up  of  a  large  number  of  program  tools  for 
interactive  analysis,  design  and  simulation  of  a  broad  class 
of  linear  control  problems.  with  LINCON,  users  can 
concentrate  on  their  specialized  applications  rather  than 
system  design  and  routine  program  development,  thereby 
saving  valuable  time. 

It  is  assumed  that  the  reader  is  familiar  with  the  basic 
concepts  of  linear  control  theory  as  may  be  obtained  from 
any  one  of  a  number  of  available  textbooks  (see  the 
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bibliography)*  fts  such,  ths  LINCON  subprograms  are 
presented  in  a  usec-aciented  fashion.  First,  their  purpose 
and  some  general  rules  that  may  apply  are  given;  then  the 
input  reguirements  are  presentei  and  the  expected  outputs 
are  described.  Examples  for  earh  are  worlcsi  out  including  a 
copy  of  the  terminal  session  and  the  final  results. 


ri 


siSSSa 


During  LINCON's  latest  development,  underlying 
guidelines  called  for  concepts  wdich  accounted  for  the  needs 
of  the  computer  system,  the  programmer,  and  the  user. 

The  guidelines  followed  luring  this  latest  development 
stage  were: 

(1)  Operation  of  the  system  should  be  in  an  on-line 

interactive  mode  such  that  data  can  be  easily  input 
to  (or  output  from)  the  system  and  readily  accessed 
for  verification,  examination,  and  processing. 

(2)  Program  development  should  be  in  a  high-level 

language  in  order  to  facilitate  software 

implementation  and  promote  machine  independence. 

(3)  The  software  should  be  aoluiar  in  structure  so  that 

programs  can  be  modified  or  inserted  without 

affecting  existing  programs. 

(U)  Programs  ehoili  be  invoiced  by  means  of  logical 

procedures  or  ccamanls  which  minimize  ir.” eraction 
time  and  which  are  user- oriented  so  that  people  can 
operate  the  system  witnout  first  becoming  computer 
experts. 
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A.  INTERACTIVE  OPERAIION 


The  applicatioas  of  LIHCDK  ace  focused  on  ii:tecacti7e 
processing.  Experiance  has  shown  that  interactive  on-line 
coDDunication  has  nan/  advantages  in  a  research  environment 
because  it  offers  rha  opportunity  to  make  observations  and 
select  alternate  courses  of  action  in  a  more  flexible  manner 
than  with  batch  procassing  [  3].  LINCON  is  organized  around 
a  collection  of  intar-relatad  command  programs,  each  of 
which  performs  a  spacifiad  function  and  can  be  executed  by 
means  of  a  simple  keyboard  initiation  sequence. 

B.  HIGH  LEVEL  LANGUAGE 

An  important  feature  in  the  design  of  LINCON  is  than  it 
was  implemented  in  a  high  level  language.  Program 
development  in  assemoly  language  is  more  time  consuming  and 
results  in  system  dependent  software. 

LINCON  is  programmed  in  AN3I  standard  FORTRAN  and 
follows  the  conventions  of  FORTRAN  IV.  FORTRAN  has  been 
found  to  be  a  useful  Language  for  several  reasons: 


(1)  since  some 

form  of 

FORTRAN  is  a 

vailaole  on  most 

com  puters. 

LINCON 

is  highly  po 

rtabie  from  one 

computer  t 

o  another. 

A  FORTRAN 

based  system  is 

helpful  for 

ia  porting 

programs  as 

well  as  exporting 
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them.  Of  coarse,  FORTRAN  compilers  don’t  all  follow 
the  sane  standards  so  tnere  can  still  be 

difficulties. 

(2)  FORTRAN  is  a  simple  aaoagh  language  that  relatively 
complex  programs  can  be  implemented  in  a  short 
period  of  time.  Most  scientific  and  research 
personnel  know  FORTRAN  sufficiently  well  to  write 
their  own  programs  if  nscessary. 

(3)  Algorithms  can  be  testai  and  implemented  in  FORTRAN 

and  later  converted  to  assembly  language  versions  if 
more  speed  and  efficisncy  are  necessary.  This 

procedure  has  the  farther  benefit  of  aiding 
portability  such  that  even  if  parts  have  been 
converted  to  assembler,  equivalent  FORTRAN  versions 
are  available. 

C.  MODOLAR  SOFTWARE 

Overall  system  flexibiity  is  achieved  by  means  of 

modularity.  LINCON  is  actually  made  up  of  a  large  number  of 

independent  command  programs.  Each  ccmmand  program  stands 
or.  its  own  with  the  ability  to  take  some  form  of  input, 
possibly  supplied  by  some  previous  command,  and  generate 
seme  form  of  output,  possibly  to  be  used  by  a  follow-up 
ccmmand. 
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D.  0SE8-0PIENTED  3PERATI3N 

An  iaportant  aspect  in  this  aodification  of  LINCON  was 
to  make  the  ccmaiaads  user-oriented  so  that  operating  the 
system  does  not  require  an  engineering  or  computer 
background.  This  was  achievai  with  a  standard  terminal 
keyboard  by  using  phonetic  characters  which  relate  to  the 
function  which  the  command  is  to  perform.  A  good 
combination  of  brerity  and  clarity  is  built  into  LINCON  to 
avoid  having  to  push  extra  buttons  on  the  keyboard  while  at 
the  same  time  preventing  ambiguity. 

S.  IMPOT  RESTRICTIONS  AND  LIMITATIONS 

Although  the  input  requirements  are  fully  described  in 
the  presentation  of  each  program,  there  are  several  input 
format  similarities  used  by  all  of  them.  For  ease  of  use, 
and,  more  honestly,  for  ease  of  programming,  most  of  the 
data  input  is  grouped  in  the  sane  arrangement. 

The  first  input  of  every  program  is  used  to  identify  the 
problem  for  reference  and  for  output  data.  A  maximum  of 
twenty  alpha-numeric  characters  can  be  used.  This 
restriction  was  not  a  system  limitation  but  a  programmer 
decision. 
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The  next  input  caamon  to  all  programs  is  the  dimension 
of  the  plant  matrix  or  &  matrix.  The  format  is  II  which 
would  normally  restrict  the  user  to  a  maximum  matrix  size  of 
9x9«  however,  due  to,  again,  a  programmer  decision,  a 
dimension  size  not  to  exceed  8  is  requested.  The  reasoning 
behind  this  was  due,  in  part,  to  the  printer.  The  NPS 
printer  is  capable  of  printing  133  characters  on  a  line. 
Since  the  output  format  to  the  printer  is  8E16.6  this 
naturally  limits  one  to  8  numbers  per  line.  Six  places  are 
normally  considered  necessary  for  good  accuracy.  A  solution 
that  would  have  lead  to  an  unlimited  matrix  size  would  have 
beer,  to  incorporate  a  "wrap  around"  routine  within  the 
program.  After  attempting  this  it  was  decided  the  results 
were  just  too  difficult  to  read. 

Matrices  are  antared  one  alament  at  a  time  beginning 
with  element  1,1  and  continui.tg  across  the  row.  The  next  row 
is  then  entered,  and  the  process  continues  until  all 
elements  have  been  entered.  After  the  matrix  is  entered,  the 
complete  matrix  is  automatically  brought  to  the  screen  for 
review  and  possible  correction.  If  a  change  to  the  matrix  is 
desired  the  user  simply  enters  the  row  number  and  column 


number  without  a  separating  conma 


For  example 


35  would 


indicate  the  eleaant  in  row  3  and  coluan  5.  After  being 
prompted  the  change  is  entered.  A  review  of  the  matrix  is 
again  screened.  The  user  is  again  prompted  for  any  possible 
changes.  This  procedure  contiiues  until  all  changes  have 
been  made. 

Any  special  reguir ements  jc  limitatious  will  be  brought 
to  the  user's  attentio.n  within  each  prograi  presentation. 
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III.  aAlfili  lEAfiSfSE  FOKlllos  and  Jias 

asspDNSE  ES2SS&as 


A.  INTHODaCTION 


In  this  chapter  four  prograis  are  discussed  which  aay  be 
used  for  the  analysis  and  design  of  linear  control  systems 
represented  in  state  variable  fora  as 


x<t)  *  45(c)  ♦  buft)  (J.  A-1) 

u  (t)  *  K[r(t)  -  jc^  j(t)  ]  (3.  A-2) 

Y  (t)  =  c'^5(t)  (3.A-3) 


The  first,  BASilAT,  is  toe  Basic  Matrix  manipulation 
program  which  is  ased  for  computing  the  determinant, 
inverse,  characteristic  polynonial,  and  eigenvalues  for  a 
square  matrix  A.  Ii  addition,  BASMAT  will  calculate  the 
state  transition  matrix  and  the  PHI (s)  matrix.  The  second 
program,  PRFIXP,  calculates  the  partial  fraction  expansion 
of  a  polynomial.  The  third  program,  ROOTS,  calculates  the 
roots  of  a  polynomial.  The  fourth  program  is  used  for 
determining  the  time  response  of  linear  control  systems. 
RTPZSP  will  determine  the  ratioaal  time  response  of  a  system 
in  closed-form  provided  that  the  input  function  r(t)  has  a 
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rational  time  response  and  that  there  be  no  repeated 
eigenvalues  in  the  coabination  of  the  systen  and  input.  It 
should  be  noted  here,  and  Kill  be  again  in  the  actual 
discussion  cf  the  program,  that  by  setting  r (t)  and  K  equal 
to  zero,  unforced  aid  open-loop  systems  may  be  studied, 
respectfully. 
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B.  B&SIC  MATRIX  PR03RAM  (BASMAT) 


Given  the  plant  matrix  A,  BASMAT  can  compute  the 


following: 


(1)  the  determinant  oE  A,  det  A 


(2)  the  inverse  of  A,  A-» 


(3)  the  characteristic  polynomial,  det(sI-A) 

(4)  the  eigenvalues  of  the  characteristic  polynomial, X 

(5)  the  state  transition  matrix,  J(t)  *  exp  (At) 

(6)  the  ?HI(s)  matrix,  f(s»  =  (sI-A)-* 


1.  Terminal  S^sei  on  Sxamplg 


This  section  contains  a  terminal  session  for  a 


specific  problem.  Commands  entered  by  the  user  are 


lower  case.  All  of  BASMAT's  capabilities  will  be  utilized 


beginning  wirh  following  plant  aarrix: 


1.0  0.  0  0.0 

A  =  0.0  -2.0  1.0 

0.0  -0.5  1.0 

(.  I 


lincon 

SXECOTION  BEGINS... 

LINCON  CONSISTS  OF  THE  FOLLOWING  S08PR0GEAMS; 

BASIC  MATRIX  HANtPOLAriON  -  <3ASMAT> 

RATIONAL  TIME  RESPONSE  -  <RrRSS?> 

STATE  VARIABLE  FEEDBACK  -  <3rVAR> 

CONTROL  LABILIT?  AND  OBSERVABILITY  -  <OBSCON> 
L'JENBERGER  OBSERVER  -  <LaEN> 

OPTIMAL  CONTROL/KALMAN  FILTERS  -  <RICAri> 

DISCRETE  TIME  SALMAN  FILTER  -  <KALMAN> 

OPTIMAL  CONTROL  -  <OPrC0N> 

PARTIAL  FRACTION  EXPANSION  -  <?RFEXP> 

ROOTS  0?  A  POLYND1IAL  -  <RD0rS> 

TO  DSE  ONE  OF  THE  SOB  PROGRAMS  ENTER  THE  NAME  BETWEEN  THE 
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SYHBOLS 

basaat 


<  >. 


BASMAT  PROVIDES  SATBIX  MANIPOLATION  TO  SOLVE 
FOR  DBTERMINAHTS,  INVERSES,  STATE  TRANSITION  AND 
PHIJS)  MATRICES,  EIGENVAL0E3,  AND  CHARACTERISTIC 
POLYNOMIALS. 

FIRST  ENTER  THE  PROBLEM  IDENTIFICATION 
(♦NOT  TO  EXCEED  23  CHARACTERS*) 
thesis  exasple 

NON  ENTER  THE  DIMENSION 


3 


THE 

1.0 

ELEMENT 

A(1,1l  = 

THE 

0.0 

ELEMENT 

A(1,2)  = 

THE 

0.0 

ELEMENT 

A(1,3|  = 

ELEMENT 

A(2,1)  = 

-2.3“' 

ELEMENT 

A(2,2|  = 

THE 

1.0 

ELEMENT 

A(2,3)  = 

0.0 

ELEMENT 

A(3,1>  = 

THE 

«0.5 

ELEMENT 

A(3,2>  =* 

THE 

1.0 

ELEMENT 

A(3,3l  * 

THE 

1  A  MATRIX 

1.COE  +  00  0.3 

0.0  -2.00E+00 

0.0  «5.30E-01 

DO  YOU  WANT  TO  CHANGE  ANT 
n 


OF  THE  A  MATRIX  (OP  TO  8) 


3.3 

1.DOE+00 
1. 30E+00 

LEMENT  OP  THE  MATRIX? 


DO  TOO  WANT  TO  CALCDLATE  THE  DETERMINANT? 

y 

THE  DETERMINANT  OF  THE  MATRIX 
-1 . 50E  +  00 

ARE  TOO  SATISFIED  WITH  THE  RESULTS? 

y 

DO  TOO  WANT  THE  DETERMINANT  PRINTED? 

y 

DO  TOO  WANT  TO  CALCULATE  THE  INVERSE? 

y 
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THE  INVERSE  OP  PHE  MATHIX 


1.00E+00  0.3  0.3 

0.0  -6.57B-31  5.57E-01 

0.0  -3,33E-01  1.33E+00 

ARE  TOO  SATISFIED  WITH  THE  RBSOLTS? 

y 

DO  TOO  WANT  THE  INVERSE  PRINTED? 

y 

DO  TOO  «ANT  TO  CALCULATE  THE  PHI (S)  MATRIX? 

y 

******** innmm*********  **********************************  *** 


y 

y 

y 


THE  MATRIX  COEPPtCIENTS  3F  THE  NOMSRATDR  OF  THE 
PHl(S)  MATRIX 

THE  MATRIX  COEFFICIENT  OF  3**2 


1.  OOE^OO 
0.  0 
0.  0 


0.3 

1.30E  +  00 
0.  3 


0.  3 

0.3 

1. 30E+00 


THE  MATRIX  COEFFICIENT  OP  3**1 


1.  OOE  +  00 

0.  0 


0.  3 
2.0 


OE+00 


0  3 

tl 30E*00 
1.30E+00 


0.0  -5.30E-01 

THE  -MATRIX  COEFFICIENT  OF  3**0 


-1.  50E+00 
0.  0 
0.  0 


I^OS+OO  -lI^OE* 


00 


5.30E-01  -2.00E  +  00 


ARE  TOO  SATISFIED  WITH  THE  RESULTS? 

A  NO  RESPONSE  WILL  GIVE  TOO  THE  OPTION  TO  MAKE 
CHANGES  TO  THE  A  MATRIX. 


DO  YOU  WANT  A  PRINTOUT  OF  THE  RESULTS? 

DO  YOU  WANT  TO  CALCULATE  THE  CHARACTERISTIC  POLYNOMIAL? 


**********************  ********************************** *** 


THE  CHARACTERISTIC  POLYNOHI AL-IN  ASCENDING  POWERS  0?  S 
1.50E  +  00  -2.  5  0E+00  3.3  1.00E+00 

ARE  YOU  SATISFIED  WITH  THE  RESULTS? 


DO  YOU  WANT  THE  CHARACTERISTIC  POLYNOMIAL  PRINTED? 

y 

DO  YOU  WANT  TO  CALCULATE  THE  EIGENVALUES? 

*i******  **  *********  ***  ************  ****  ****************** 


THE  BIGENVALOES  3P  THE  A  MATRIX 
REAL  PART  IBAS.  PART 
8.23E-01  0,3 

-1.82E  +  00  0.3 

1.00E  +  00  0,3 

ARE  YOO  SATISFIED  WITH  THE  RESOLTS? 

y 

DO  YOD  WANT  THE  El  GEM? ALOES  PRINTED? 

y 

DO  YOO  WANT  TO  CALCOLATE  THE  STATE  TRANSITION  MATRIX? 

y 

4i3k*4ii|t4i4>4i4t*4>4>4>*4>4t4t4>4<4>#4i  ***■*******■¥******  *******  *** 

THE  ELEMENTS  DP  THE  STATE  TRANSITION  MATRIX 

THE  MATRIX  COEFFICIENT  OP  EXP  (  8.23E-01)T 

0.  0  0.3  3.3 

0.0  -6.59E-32  3.  7BE-01 

0.0  -1.39B-01  1.37S4-00 

THE  MATRIX  COEFFICIENT  OF  EXP(  -1.82E«-00)T 


0.  0 
0.  0 
0.0 


iISte+OO  -3l?SE>01 
1.39E-01  -6,59E-02 


y 

y 


THE  MATRIX  COEFFICIENT  OF  EXP  (  1.00E«-30)T 

1.g0E  +  00  .JlJsE-Oe  ?^?6B-06 
0.0  -5.36E-07  1.37E-06 

ARE  YOU  SATISFIED  WITH  THE  RESOLTS? 

A  NO  RESPONSE  MILL  GI? E  THE  DPTION  TO  MAKS 
CHANGES  TO  THE  A  MATRIX. 

DO  YOU  WANT  A  PRINTOUT  OF  THE  RESULTS? 

THIS  CONCLUDES  THE  BASIC  MATRIX  MANIPULATION  PROGRAM 
(BASMAT) . 

ANALYSIS  IS  COMPLETE.  DO  TOO  WANT  TO  RUN  LINCON  AGAIN? 
n 

LINCON  IS  NOW  TERMINATED. 

The  computer  results  are  shown  in  Appendix  A. 
Interpretation  of  the  determinant,  inverse  and  eigenvalues 
are  straightforward.  The  PHI(st  matrix  is  a  combination  of 
the  numerator  coefficients  and  characteristic  polynomial. 
The  results  can  be  interpreted  as 
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and,  sialarily,  tha  fifth  tarm  of  tha  state  transition 
matrix  is 

f(t)  =  -0.067exp(D-82t)  ♦1.07axp(-1.82t> -2.26x10-‘9xp  (t) 
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C.  PARTIAL  FRACTION  EXPANSION  PROGRAM  (PRFEXP) 

PHFEXP  car  calculate  a  partial  fractian  expansion  given 
a  rational  ratio  of  tvo  polynoaials  in  the  fora 

N  (e) 

G(3)  =  K - -  (3.C-1) 

0(3) 

atiere  K  =  the  input  function  gain, 
a  constant 

N  {s(  =  n,  +  n^s  ♦  nj  s2  s' 

and  D(s)  =  d,  ♦  s  +  dj  s2  ♦  ...  + 

The  nunerator  and  denominator  coefficients  must  be  arranged 
so  that  the  coefficiants  of  s  and  s  are  each  unity  with 
p  >  q  >  0.  The  outputs  of  the  program  are 

(1)  the  numerator  gain 

(2)  the  numerator  polynomial  and  its  roots;  roots  are 
considered  equal  if  th=ir  real  and  imaginary  parts 
do  not  differ  by  more  tian  0.005. 

(3)  the  denominator  polynomial  and  its  roots;  multiple 
roots  are  listed  once,  along  with  its  multiplicity 

(U)  the  numerator  ccefficients  (residue  matrix)  are 
listed  in  the  same  order  as  the  denominator  roots; 
the  first  coefficient  in  any  row  of  the  matrix  is 
the  first-order  term  coefficient. 


2'4 


1.  iassion 


The  partial  fraction  axpansion  of  the  following 
rational  polynomial  is  to  be  performed: 


G(S) 


14  (s2+2)  (s+1> 
2s^  +  S  s2+2s*5 


Putting  the  polynomial  in  a  usaole  form  yields 


G  (S) 


7  (ss*  s2+2s  +  2> 
+  2. 5 


lincon 

EXECUTION  BEGINS... 

LINCON  CONSISTS  OF  TNE  FOLLOWING  SUBPEOGHAMS: 

BASIC  JIATRIX  MAKIPOLAI  ION  -  <BAStlAT> 

NATIONAL  TIME  RESPONSE  -  <RrRESP> 

STATS  VARIABLE  FEEDBACK  -  <3rVAR> 

CONTROL LABILITX  AND  OBSERVABILITY  -  <OBSCON> 
LUENBEHGEH  OBSERVER  -  <LUBN> 

OPTIMAL  CONTROL/KALMAN  FILTERS  -  <EICAri> 

DISCRETE  TIME  KALMAN  FILTER  -  <KA1MAN> 

OPTIvMAL  CONTROL  -  <OprCON> 

PARTIAL  FRACTION  EXPANSION  -  <PRFEXP> 

ROOTS  OF  A  POLYNOMIAL  -  <ROOrS> 

TO  USE  ONE  OP  THE  SUBPROGRAMS  ENTER  THE  NAME  BETWEEN  THE 
SYMBOLS  <  >. 
prf exp 

PREEX?  IS  USED  TO  DETERMINE  THE  PARTIAL 

FRACTION  EXPANSION  OF  THE  RATIO  OF  TWO  POLYNOMIALS. 

FIRST  ENTER  THE  PROBLEM  IDENTIFICATION 
(♦NOT  TO  EXCEED  20  CHARACTERS*)  . 
rhesis  exaipie 

ENTER  THE  INPUT  FUNCTION  3AIS--K. 


THE  JUNCTION  GAIN 

DO  YOU  WANT  TO  CHANGE  THE  VALUE  OF  THE  GAIN? 
n 

ENTER  THE  NUMERAIOR  BY  POLYNOMIAL  COEFFICIENT  OR 
FACTORED  ROOT  FORM.  FIRST  ENTER  EITHER  THE  LETTER  ? 
FOR  POLYNOMIAL  COEFFICIENT  FORM  OR  THE  LETTER  ?  FOR 
FACTORED  ROOT  FORM . 


ENTER  THE  NOMERArDR  POLYNOMIAL  ORDER. 

3 

THE  POLYNOMIAL  COEFFICIENTS  .1UST  BE  ENTERED  IN  ASCENDING 
ORDER  OF  S. 

♦♦WARNING — THE  HIGHEST  ORDE3  COEFFICIENT  MOST  BE  ONITY.^ 
DO  YOO  NEED  TO  CHANGE  THE  INPOT  FONCTION  GAIN  TO  SATISFY 
THIS  REQOIREMENT? 


ENTER  THE  POLYNOMIAL  COEFFICIENTS  IN  ASCENDING 
ORDER  OF  S. 

COEFF( 1)* 


COEFP(2)* 


COEFF(3)* 


COEFP(4)* 


NOMEPATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S 
2.COE+00  2.00E+00  1.00E  +  00  l.OOE+00 


p 

u 


NOMEPATOR  ROOTS  ARE 

REAL  PART  IMAG  .  PART 
0.0  -1.41E+00 

0.0  1.41E+00 

-1.C0E  +  00  0.3 

ENTER  THE  DENOMINATOR  BY  POLYNOMIAL  COEFFICIENT  OR 
FACTORED  ROOT  FORM.  FIRST  ENTER  EITHER  THE  LETTER  ?  FOR 
POLYNOMIAL  COEFFICIENT  FORM  OR  THE  LETTER  P  FOR  FACTORED 
HOOT  FORM. 


ENTER  THE  DENOMINATOR  POLYNOMIAL  ORDER. 


THE  POLYNOMIAL  COEFFICIENTS  MUST  BE  ENTERED  IN  ASCENDING 
ORDER  OF  S, 

♦•HARNI NG--THE  HIGHEST  ORDER  COEFFICIENT  MUST  BE  jNITY.* 
DO  YOU  NEED  TO  CHANGE  THE  INPUT  FUNCTIDN  GAIN  TO  SATISFY 
THIS  REQUIREMENT? 


ENTER  THE  POLYNOMIAL  COEFFICIENTS  IN  ASCENDING 
ORDER  OF  S. 

COEFF(1  )* 


COEFF(2  )* 


COEFF(3  )* 

0 


C0EF7  (4  )  = 
COEFF  (5  )  = 
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DENOMINArOR  CDEPPICIENTS  -  IN  ASCENDIN3  POWERS  OF  S 
2.50E+00  1.00E+00  3  .0  3.00E+00  1.00E+00 


DENOMINATOR  ROOTS  ARE 
REAL  PART  IMA3.  PART 
-3.019E  +  00  0.3 
-9.  128E-01  0.0 
4.658E-01-8.308E-3  1 
4.658E-01  8.308E-31 


MULTIPLICITY 

1 

1 

1 

1 


RESI^UE^MAJglX  -  REAL  PART 

3!  17E-01 
4.  36E-01 
4.  36E-01 


iha: 


PART 


RESIDU^MAJ^IX 

Ol  0 

2.  C6E+00 
-2.  06E+00 

THIS  CONCLUDES  THE  PARTIAL  FRACTION  EXPANSION  PROGRAM 
(PRFEXP) . 


DO  YOU  WANT  TO  RUN  THE  PROGRAM  AGAIN? 


ANALYSIS  IS  COMPLETE.  DO  YOU  WANT  TO  RUN  LINCON  AGAIN? 
n 

LINCON  IS  NOW  TERMINATED. 

The  resul'^s  are  shown  in  Appaaiix  B.  Interpretation  of 
these  resulxs  are; 


_  s-5.81*j0  3-0.32  3-0.44-j2.06  s-0.a4*j2.06 

“  "s+iroi  *  s+o.ii’  ’’  s-oTaYT jo.83”  *  s-oraT-jo.ai 
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D.  POLYNOMIAL  ROOIS  PROGRAM  (ROOTS) 

ROOTS  can  calculate  the  roots  of  a  polynoaial  of  degree 
less  than  or  equal  to  eight.  Gi/ea  a  pclynoaial  in  the  fora 

P(s)  =  p.  ♦  PiS  «•  Pi  s2  V  ...  ♦  s'  (3.D-1) 

the  coefficient  of  3^  must  be  unity.  The  program  output 
lists  the  polynomial  coefficiants  for  reference  and  the  real 
and  imaginary  roots. 

1.  lerjdiial  Session  Example 

The  following  polynomial  is  to  be  factored: 

P(s)  =  “5s3  15s2  -  25. Ss"  4 

Putting  the  polynomial  in  the  required  fora  yields 

?  (s)  =  S3  -  3s2  «•  5. 13  -  0.8 


iir.con 

EXECUTION  BEGINS _ 

LINCON  CONSISTS  OF  IMS  FOLLOiJINS  SUBPROGRAMS; 

BASIC  MATRIX  MANIPULATION  -  <5ASMAT> 

RATIONAL  TIME  RESPONSE  -  <RrRESP> 

STATE  VARIABLE  FEEDBACK  -  <3rVAR> 

CONTROLLABILITY  AND  OBSERVABILITY  -  <OBSCON> 
LUENBERGER  OBSERVER  -  <LUEN> 

OPTIMAL  CONTROL/KALMAN  FILTERS  “  <RICAri> 

DISCRETE  TIME  KALMAN  FILTER  -  <KALMAN> 

OPTIMAL  CONTROL  -  <OPrCON> 

PARTIAL  FRACTION  EXPANSION  -  <PRFEXP> 

ROOTS  OF  A  POLYNOMIAL  -  <ROOrS> 

TO  OSS  ONE  OF  THE  SUBPROGRAMS  ENTER  THE  NAME  BETWEEN  THE 

SYMBOLS  <  >. 

rccts 
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5S851e"MFIhIS  5i''ijS2E 

FIHST  ENTER  THE  PROBLEM  IDENTIFICATION 
(NOT  TO  EXCEED  2D  CHARACTERS)  . 
thesis  eza  ople 

ENTER  THE  ORDER  DP  THE  SYSTEM  (OP  TO  8). 


ROOTS  OF  A  POLYNOMIAL 

ENTER  THE  POLYNOMIAL  COEFFIZIENTS  IN  ASCENDING 
ORDER  OF  S. 

♦♦WARNING  — THE  HIGHEST  ORDER  COEFFICIENT  MOST  BE  ONITY.^ 
COEFP(1)= 


COEFF(2)* 

5.1 

COEPP(3)= 


COEFF(U)* 


POLYNOMIAL  COEFFICIENTS  -  IS  ASCENDING  POWERS  OF  S 
8.00E-01  5.10E4-00  3.DDE  +  00  1.0DE  +  00 

THE  ROOTS  ARE  REAL  P\Rr  IMAG.  PART 

-1.41E*-00  1.62S+3D 

"1.aiE4-00  -1.62E4-00 
-1.74S-01  0.0 

THIS  CONCLODES  THE  POLYNOMIAL  ROOTS  PROGRAM. 

DO  YOO  WANT  TO  RON  THE  PROGRAM  AGAIN? 


ANALYSIS  IS  COMPLETE.  DO  YOO  WANT  TO  RON  LIKCON  AGAIN? 


LINCON  IS  NOW  TERMINATED. 

The  results  shown  i?.  Appendix  C  are  sasily 


inter pretad . 
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E.  RATIONAL  TIME  RESPONSE  PROGRAM  (RTRESP) 

The  time  response  in  closed  form  of  a  linear  control 
systea  described  by  Eqs.  (3.A-1),  (3.A-2)  and  (3.A-3)  is 

calculated  by  this  program.  The  user  must  define  the 
initial  conditions  5(3)  and  tie  rational  Laplace  transform 
of  the  scalar  forring  function  r(t).  The  theoretical 
concepts  involved  in  the  development  of  the  computer  codes 
are  described  by  Melsa  C  ^  1* 

Inp ut  Restrict  ions  And  Limitations 

Enter  the  following  ia  F  format  (floating  point)  : 

1)  the  elements  of  the  plant  matrix  A 

2)  the  control  vector  b 

3)  the  output  vector  c 

4)  the  feed  bacicback  coefficent  vector  jc 

5)  the  controller  gain  < 

and 

6)  the  initial  conditions  vector  x(0). 
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-J-MJ.  iJ  .JJHU 


The  rational  Laplace  transform  of  the  input  function 
must  be  in  the  form 

'f  M(S» 

oL.  [r(t)  ]  *  8(S»  »  G -  (3.  E-1) 

D(3» 

i/here  G  »  the  input  function  gain, 
a  constant 

N  (s)  s  n,  ♦  n^s  «•  n^sz  ♦  ...  +  s*" 
and  D(s)  =  d,  djS*'  djsz  ♦  ...  ♦ 

The  numerator  and  denooinaxor  coefficients  must  be  arranged 
so  that  the  coefficients  of  s'  and  s^  are  each  unity  with 
p  >  g  >  0.  Due  to  programmin j  limitations  it  is  necessary 
that  the  order  of  the  system  plus  the  dimension  of  D (s)  be 
less  than  or  equal  to  ten. 

Upon  entering  the  input  function  gain,  the  user  next 
has  the  option  to  enter  the  numerator  and  denominator  in 
either  polynomial  coefficient  fora  or  factored  root  form. 
With  the  factored  root  form  entsc  the  real  part  of  the  root 
as  negative  if  it  lies  in  the  Left  half  plane  and  just  the 
magnitude  of  the  imaginary  part. 

Session  3xamBl= 

This  section  contains  a  terminal  session  for  a 
specific  example.  Coaaands  entered  by  the  user  are  in  lower 
case.  The  RTR3SP  program  will  determine  the  time  response 
form  of  the  closed-loop  system 
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I(t) 


'  1.0 
=  0.0 
0.0 
L 


0.0 

-2.3 

-0.5 


0.0' 

1.0  x(t) 

1.0 

i 


0.0 

♦  0.0  u(t| 
1.0 
L  J 


u(t)  =  3.2  {r(t)  - 
I(t)  =  [1.0  0.4 

if  the  Laplace  transf 

a(s)  =  0.5  - 


[1.0  1.0 
1.0]  5(t) 
arm  of  the 

3  ♦  1.0 

s*~-~2ro 


0.0]  x(t)» 


input  function  is 


and  the  initial  conditions  ars  zero. 

Note  here  that  the  open-loop  rational  time  response 
may  be  calculated  by  setting  the  feedback  coefficient  vector 
k  equal  to  a  zero  vector.  Also,  if  only  an  initial 
condition  response  is  desired,  the  input  function  gain  G  is 
set  to  zero. 


linccn 

EXECUTION  BEGINS... 

LINCON  CONSISTS  OF  TSE  FOLLOWING  SUBPROGRAMS: 

BASIC  MATRIX  MANIPULAIION  -  <3ASMAT> 

RATIONAL  TIME  RESPONSE  -  <arRESP> 

STATS  VARIABLE  FEEDBACK  -  <5TVAS> 

CONTROLLABILITY  AND  OBSERVABILITY  -  <OB5CON> 
LUZN3ERGER  OBSERVER  -  <LUSN> 

OPTIMAL  CONTROL/KALMAN  FILTERS  -  <EIC^.ri> 

DISCRETE  TIME  KALMAN  FILTER  -  <KALMAN> 

OPTIMAL  CONTROL  -  <CPrC0N> 

PARTIAL  FRACTION  EXPANSION  -  <PEPEXP> 

ROOTS  OF  A  POLYNOMIAL  -  <RO0TS> 

TO  USE  ONE  OF  THE  SUBPROGRAMS  ENTER  THE  NAME  BETWEEN  THE 

SYMBOLS  <  >. 

rtresp 
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T 


1- 

0. 

1. 


1 . 
1 , 

0. 


n 

3. 


0. 

0. 

G .  0 


ENTER  THE  ELEBENIS  OF  THE  DOTPOT  VECTOR--C. 

jC(1)» 

^C(2)» 

jC(3)» 

THE  C  MATRIX  (OOTPUT  VECTOR) 

1. QQE+00 
4,  COE-01 
1 .  OOE  +  00 

DO  YOU  WANT  TO  CHANGE  ANY  ELEMENT  (S)  OF  THE  MATRIX? 

ENTER  THE  ELEMENTS  OF  THE  FEEDBACK  COEFFICIENT  VECTOR 
FDBG. 

^FDBG  (1)  = 

jFDBG  (2)  * 

^FDBG  (3)  = 

THE  FEEDBACK  COEFFICIENT  VECTOR 

1. OOEfOO 
1.  COE  +  00 
0.0 

DO  TOO  WANT  TO  CHANGE  ANY  ELEMENT  (S)  OF  THE  VECTOR? 

^ ENTER  THE  CONTROLLER  GAIN--K. 

THE  CONTROLLER  GAIN 
3. 20E>00 

DO  YOU  WANT  TO  CHANGE  THE  VALUE  OF  THE  SAIN? 

ENTER  THE  ELEMENTS  OF  THE  INITIAL  CONDITIONS  VECTOR  - 
X  (0)  . 

^XO(1)  = 

jXO(2)* 

XO(3)* 

INITIAL  CONDITIONS  VECTOR-X (0) 

0.  0 
0.  0 
0.  0 
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DO  TOO  WiST  TO  CHMIGE  ANY  ELEMENT  (S)  OF  THE  VECTOR? 


ENTER  THE  GAIN  FDR  THE  RATIONAL  LAPLACE 
TRANSFORM  OF  THE  INPOT  FONCTION. 


THE  INPOT  FONCriDN  GAIN  = 

5.  OOE-01 

DO  TOO  WANT  TO  CHANGE  THE  VALOE  OF  THE  GAIN? 


ENTER  THE  NOMERATDR  BY  POLYNDHIAL 
COEFFICIENT  OR  FACTORED  ROOT  FORM. 

FIRST  ENTER  EITHER  THE  LETTER  P  FOR 
POLYNOMIAL  COEFFICIENT  FORM  DR  THE  LETTER  ?  POE 
FACTORED  ROOT  FORM. 


ENTER  THE  NUMERATOR  POLYNOMIAL  ORDER. 


ENTER  THE  REAL  PART  OP  THE  RDDT. 

-1.0 

ENTER  THE  MAGNITUDE  OP  THE  IMAGINARY  ROOT. 

0.0 

NUMERATOR  POLINDMIAL  OF  R (S> -ASCENDING  POWERS  OF  S 

1.00E  +  00  1.00E+30 

NOMEBATOR  ROOTS  ARE 
REAL  PART  IMA3.  PART 
-I.OOE-t-00  O.D 

ENTER  THE  DENOMINATOR  BY  POLYNOMIAL 
COEFFICIENT  OR  FACTORED  ROOT  FORM. 

FIRST  ENTER  EITHER  THE  LETTER  P  FOB 
POLYNOMIAL  COEFFICIENT  FORM  DR  THE  LETTER  ?  FOR 
FACTORED  ROOT  FORM . 

P 

ENTER  THE  DENOMINATOR  POLYNOMIAL  ORDER. 

2 

THE  POLYNOMIAL  COEFFICIENTS  MUST  BE 
ENTERED  IN  ASCENDING  ORDER  OP  S. 

WARNING~THE  HIGHEST  ORDER  COEFFICIENT  MUST  BE  UNITY 
DO  YOU  NEED  TO  CHANGE  THE  INPUT 
FUNCTION  GAIN  TO  SATISFY  THIS  REQUIREMENT? 
n 

ENTER  THE  POLYNOMIAL  COEFFICIENTS  IN 
ASCENDING  ORDER  OP  S. 

^  CO(1)* 

-2.0 

.  CO(2)  = 


.CO  (3)  = 
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DEHOHIN&TOR  POLINOMIAL  OP  5(5)  -  ASCENDING  POWERS  OP  S 


-2.00E  +  00  0.0 


1.00E+00 


DENOMINATOR  ROOTS  ARE 
REAL  PART  IMA3.  PART 
1.41E  +  00  0.0 

-1.41E+00  0.0 

*jHt*****0*******^**iHt*  *4i4i4i 41  ««** ******* 4>*«***4i#***4i  *** 

THE  TIME  RESPONSE  OF  THE  STATE  X  (T) 

THE  VECTOR  COEFFICIENT  OF  EXP  (-5 .  OE-01 )  T*COS  ( 1 . 2E  +  0  0)  T 
0.0  -3.70E-01  -9.72E-01 

THE  VECTOR  COEFFICIENT  OF  EXP  (-5.  OE-01 )  T*SIN  (1 . 2E  +  0  0)  T 
0.0  -3.47E-01  -7.50E-02 

THE  VECTOR  COEFFICIENT  OF  EXP  (  1.00E«-00)T 

0.  0  0.  0  2.  38E-06 

THE  VECTOR  COEFFICIENT  OF  EXP  (  1.41E+00)T 

0.  0  2.57E-01  9.12S-01 

THE  VECTOR  COEFFICIENT  OF  EXP  (  -1-41E  +  00)T 
0.0  1.03E-01  3.00E-02 

*41  ***********  4i*4t4i*4c4i*4t4t  «*  41**  41****4141  41****  4c  «*4t4i«*  *** 

THE  TIME  RESPONSE  OF  THE  OOTPUT  Y  (T) 

THE  COEFFICIENT  OF  EXP(  -5 . 00  S-0  1 )  T*C0S  (  1.20E  +  00)T 

-1  .  12E  +  00 

THE  COEFFICIENT  OF  EX?(  -5.  OOS-0  1)  T*SIN  (  1.20E  +  00JT 

-2.  14E-01 

THE  COEFFICIENT  OF  EXP  (  I.OOE  +  OO)! 

2.  68E-06 

THE  COEFFICIENT  OF  EXP(  1.41S+00)T 
1.  02S  +  00 

THE  COEFFICIENT  OF  EXP  (  -1.4lS4-00)T 
1 .  01E-01 

***********  ***********  ****  ******************************  *** 

THIS  CONCLUDES  THE  RATIONAL  TIME  RESPONSE  PROGRAM 

(RTRE3?) 

ANALYSIS  IS  COMPLETE.  DO  YOU  WANT  TO  RUN  LINCON  A3  AIN 


LINCON  IS  NOW  TERMINATED. 


35 


The  computec  results  shown  in  Appendix  D  are 


interpreted  as; 


*1  (t) 

Xi  (t) 


=  0.0 


-0 . 37exp(' 0,5t) cos  n . 2t) -0. 35exp (-0 . 5t) sin (1 . 2t) 
♦  0. 67exph -“Itl  *0.  tOexp  ^ 


X.  (t)  =  -Q .  97expj[- 0.  5t)  cns  (1 . 2t) -0.  OOexp  (-0, 5t)  s^n  (1  '**’ 
^  ♦2.7x10-Sexp<t)*0.5lexp(l.U1t)  ♦■5.06exp(-^.«1 


y(t)  *  -1 .  12exp  (- 0. 5t)  cos  (1 . 2r) -0 . 2  Isxp  (-0.  5 1)  sin  f  1 . 2t) 
♦  2.7x13-*exp  (t{  ♦1.d2axp  (l.U1r)  <'0.  1exp(-1.4k) 
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IV.  MOD^HN  COHIROi  PROGRAMS 

A.  INTRODOCTION 

In  this  chapter  six  pcograas  are  presented  and  discussed 
which  may  be  used  Enr  the  analysis  and  design  of  control 
systems. 

The  first,  0BSC3M,  is  used  to  find  the  observability 
index  and  controllability  of  a  system.  The  next  two  programs 
are  used  to  design  Kalman  filters.  RICATI  and  KALMAN  can 
find  the  feedback  ani  control  gains  necessary  to  optimize  a 
function  for  either  continuous  or  discrete  systems, 
respectively.  The  last  three  programs  may  be  used  to  design 
optimal  linear  control  systems.  srVAR  is  particularly  useful 
in  the  design  of  linear  state  variable  feedback  control 
systems.  It  may  he  used  to  calculate  both  open-  and 
clcsed*loop  transfer  E unctions  and  also  has  the  ability  to 
desiqn  a  closea-loop  system  from  desired  transfer  function 
specifications.  L(JSN  is  used  to  design  a  combined, 
reduced-order,  obser ver- concr o II er  to  achieve  a  desired 
clcsed-loop  transfer  function  from  a  system  where  some  of 
the  states  are  inaccesible.  OPTION  will  fninimize  a  given 
cost  function  producing  a  scalar  control. 
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8.  OBSEByABILTIY  INDEX  AND  CDNT ROLLABILIIT  PROGRAM  (OBSCON) 
This  program  is  used  to  determine  the  observability 
index  and  controllability  of  the  linear  system 

x(t)  =  Ax(t)  ♦  ga(t>  (4.3-1) 

I(t)  =  cx(t)  <4.3-2) 

where 

X  =  state  vector  (n- dimensional  vector) 
u  *  control  vector  (1-dimensional  vector) 

2  =  output  vector  (m-dimeasional  vector) 

A  *  n  X  n  marrix 

B  =  n  X  1  matrix 

C  =  m  X  n  matrix 

The  observability  index  is  defined  as  the  minimuni 
in-^eger  such  that  the  matrix 

rc,4  c,... ,  (k  )''-»c] 

has  rank  n.  The  above  system  is  said  to  be  controllable  at 
a  given  initial  time  if  it  is  possible,  by  using  an 
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unconstrained  control  vector,  to  force  the  sysrem  from  an 
initial  state  of  5(0t  to  sons  other  stats  in  a  finite  time 
interval  [  4], 

If  the  user  desires  just  the  observability  index  to  be 
calculated,  enter  ths  B  oatrix  is  a  zero  matrix.  Likewise, 
if  just  the  controllability  of  ths  system  is  desired,  snter 
the  C  matrix  as  a  zero  matrix.  The  unobservaola  or 
uncontrollable  system  responsa  is  then  of  course  ignored. 

I*  ISEUt  Restrict  ions 

Other  than  the  limitations  to  the  problem 
identification,  (23  characters)  and  system  size  (8x8)  there 
are  no  restrictions  to  OBSCDN.  k  reminder  here  may  be 
helpful,  however.  Ssmember  to  enter  your  elements  of  ths 
matrices  in  F  forman. 


2-  Ter  minal  Session  Example 

This  section  contains  a  session  for  a  specific 
problem.  Commands  entered  oy  che  user  are  in  lower  case.. 
The  following  system  is  to  be  tested: 


x(t) 


0.0 

=  -1  .0 

0.0 

L 


1,0 

-0.  5 

0.0 


0.0 

1.0  x(tt 

1.0  ■ 

1 


2.0 

0.0 

0.0 


1.0 

4(i) 

0.0 

I 


n 


*^0.0 

lit)  =  1.0 

-2.0 

L 


o.o' 

0.0 

0.0 

J 


THE  OBSERffiVBrLrTY  INDEX  AND 
A  SYSTEH. 


OBLSM  IDENTIFICATION 
CHASACTEas*) 


lincon 

EXECOTION  BEGINS... 

LINCON  CONSISTS  OF  THE  FOILOHING  SDBPROGHAMS: 

BASIC  aATHIX  MANIPOLATION  -  <BASMAT> 

RATIONAL  TIME  RESPONSE  -  <HrRESP> 

STATS  VARIABLE  FEEDBACK  -  <3rVAR> 

CONTROL  UBILITY  AND  OBSERVABILITY  -  <OBSCON> 

LOENBEHGER  OBSERVER  -  <LOEN> 

OPTIMAL  CONTHOL/KALMAN  FILTERS  -  <aiCAri> 

DISCRETE  TIME  KALMAN  FILTER  -  <KALMAN> 

OPTIMAL  CONTROL  -  <OPrCOM> 

PARTIAL  FRACTION  EXPANSION  -  <PEFEXP> 

ROOTS  OF  A  POLYNOMIAL  -  <ROOPS> 

TO  OSE  ONE  OP  THE  SOBPROGRAMS  ENTER  THE  NAME  BETWEEN  THE 
SYMBOLS  <  >. 
obscon 

OBSCON  DETERMINES  THE  OBSERVABILITY  INDEX  AN 
CONTHOLUBILITY  OF  A  SYSTEM. 

FIRST  ENTER  THE  PROBLEM  IDENTIFICATION 
(♦NOT  TO  EXCEED  2D  CHARACTERS*) 
thesis  exasple 

^  NOW,  ENTER  THE  ORDER  OF  THE  SYSTEM  (UP  TO  8) 

ENTER  THE  ELEMENTS  OF  THE  PLANT  MATRIX--A. 

n 

0.0 

,  „A(1,2)* 

1.0 

„  A(1,3)  = 

0,0 

,  A  (2,1)  = 

-  1.0 

,  ,A(3,1)  = 


A  (3,2)  = 


A  (3,3)  = 


THE  i  MATRIX  /PLANT  MATRIX) 

0.0  1.00E+00  0.0 

-1.00E+00  -5.  OOE-01  1.00E+00 
0.0  0.0  1.00E+00 

DO  TOO  WANT  TO  CHANGE  ANY  ELEMENT  OP  THE  MATRIX? 

n 

ENTER  THE  NUMBER  OF  COLUMNS  OF  THE  B  MATRIX. 


ENTER  THE  ELEMENTS  OF  THE  B  MATRIX. 


jB(2,1)  = 
jB(2,2)  = 

,B(3,1)  = 
B(3,2)  = 


THE  B  MATRIX 

l:88i:SS 


DO  YOU  WANT  TO  CHANGE  ANY  ELEMENT  OP  THE  MATRIX? 


ENTER  THE  NUMBER  OF  OUTPUTS. 


ENTER  THE  ELEMENTS  OF  THE  C  MATRIX. 


I*E^^  THE 


jC(1,2)  = 

C  (1,  3)  = 


C  (2,  1) 


C(2,2)  = 


C  (2,  3) 


,  C(3,  1)  « 

■2.0 

,  C(3.2)  * 


C  (3,  3)  * 


THE  C  MATRIX 


0.0  1.00E+00  0.0 

1.00E+00  I.OOE^OO  0.0 
-2.00E>00  I.OOE+'OO  0.0 

DO  YOU  WANT  TO  CHftllGE  ANY  ELEMENT  OF  THE  MATRIX? 


OBSERVABILITY  INDEX  *  2 

THE  SYSTEM  (A,  B|  IS  UNCONTROLLABLE 

^nim*m**  *******  *******  **********************************  *** 

DO  YOU  WANT  TO  RUN  OB3COM  A3AIN? 


THIS  CONCLUDES  THE  OBSERVABILITY  INDEX/CONTSOLLABILITY 
PROGRAM  (OBSCON) 

ANALYSIS  IS  COMPLETE.  DO  YOJ  WANT  TO  RUN  LINCON  AGAIN? 


LINCON  IS  NOW  TERMINATED. 

The  complete  results  presented  in  Appendix  E  should 


be  self-explanatory  to  the  user. 


C.  OPTIMAL  COMTROL/KALMAM  FILTER  PROGRAM  (RICATI) 


The  transient  solution,  mitrix  gains,  to  the  Hiccati 
differential  equatioas  for  tha  state-regulator  controller 
and  the  continuous  Xalaan  filter  are  calculated  b7  the 
RICATI  program. 

Given  a  state- regulator  problem  with  the  linear, t ime- 
in variant  system 

i  (t)  =  Ax  (t)  ♦  33  (t)  (4.C-1) 

I(t)  =  Cx(t)  (4.C-2) 

the  Riccati  equation  is  defined  as 

£(t)  =  -E(t)A-rg(t)*g  (t)BR-*B'^P{t)-Q  (4.  C-3) 


with  P  (to )  as  the  boundary  coalition.  If  u(t)  is  not 
constrained,  a  gain  aatrix  can  be  found  such  that  the  cost 
function 

J=1/2r  (t^  )E(*-4  )5(t4t  ]♦  V2/[x'^(t>  2X(t)  tu''(t)  Ru(t)  ]dt  (U.  C-4) 

% 

is  minimized  [  5].  Such  a  gain  matrix  is  defined  as 


S<.(t)  »R-»S^P{t| 


(4.  C-5) 


and 


4  4 


a(t) 


'Sg  (t)x(tl 


(“.c-e) 


The  transient  solution  is  solve!  by  the  computer.  Note  that 
the  output  from  the  computer  for  the  gain  matrix  does  not 
include  the  negative  sign  of  the  feedback  loop. 

Given  a  continous  Kalman  filter  problem  with  the  linear, 
time-invariant  system 

kit)  =  U(t)  ♦  l!t(t)  (4.  C-7) 

2(t)  *  Cx  (t)  ♦  V<t)  (4.  C-8) 

where  w(t)  ,  the  random  process  forcing  input  and  v(t)  ,  the 
measurement  noise,  have  covariaace  matrices  of  R,  the  random 
input  "covariance  matrix,  and  2»  'the  measurement  noise 
covariance  natrix,  respectfully,  the  Riccati  equation  is 
defined  as 

f(t)  =  AP(t)+P(t)i  ♦82B^-?(t)  C^a-»C|{t)  (4.  C-9) 

with 

I(to)  =  E([x(t,)-|  ]:r(*o>-I  f\  (4.Z-10) 

as  the  initial  condition  ooundary.  The  gain  matrix  found  is 
defined  as 

G^(t)  =  I-‘Cg(t) 
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(4.:-ii) 


licjiail  ustkoa 


This  section  contains  a  session  for  the  specific 
example  presented  by  Melsa  [  1].  Commands  entered  by  the 
user  are  in  lower  case.  Given  the  second  order  linear  system 


S(t)  = 


-1.0 


’’i.o  o.o'' 

s(t»  *  ^  ^  „  “(t) 

0.0  1.0 

C  J 


I(t) 


1.0  0.0 

X  (t| 

0.0  2.0 

L  J 


determine  the  optimal  transient  response  control  and  filter 
gains  for 


2  *  2 


’’l.O  i.o’ 


1.0  1.0 

L  J 


?cr  the  control  option  an  initial  time  of  0.0  and  a  final 
time  of  10.0  is  used.  For  the  fil-sr  option  an  initial  rime 
of  0.0  and  a  final  tine  of  5-0  is  used.  Also,  for  the  filter 
option,  the  initial  condition  matrix  is  chosen  to  be 
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2(*-*) 


^0.0  Q.o' 

0,0  0.0 

L  J 


For  both  options  ten  equally  spaced  values  of  G  and  G  on 
the  time  interval  t  <  t  <  t  ars  used,  i.e.,  NPOIMT  is  set 
equal  to  10. 


lir.con 

EXECUTION  BEGINS... 

LINCON  CONSISTS  OF  THE  FOLLOWING  SUBPROGRAMS; 

BASIC  MATRIX  MANIPULATION  -  <BASMAT> 

RATIONAL  TIME  RESPONSE  -  <RrRESP> 

STATE  VARIABLE  FEEDBACK  -  <3rVAR> 

CONTROLUBILITX  AND  OBSERVABILITY  -  <03SCON> 

LUENBERGER  OBSERVER  -  <LUEN> 

OPTIMAL  CONTROL/KALMAN  FILTERS  -  <RICAri> 

DISCRETE  TINE  KALMAN  FILTER  -  <KALHAN> 

OPTIMAL  CONTROL  -  <OPrCON> 

PARTIAL  FRACTION  EXPANSION  -  <PRFEXP> 

ROOTS  OF  A  POLYNOMIAL  -  <RO0TS> 

TO  USE  ONE  OF  THE  SUBPROGRAMS  ENTER  THE  NAME  BETWEEN  THE 
SYMBOLS  <  >. 
r icati 

RICATI  DETERMINES  THE  TRANSIENT  SOLUTION  FOR  THE  RICCATI 
EQUATION.  FIRST  ENTER  THE  PROBLEM  IDENTIFICATION 
(*NOT  TO  EXCEED  20  CHARACTERS*) 
thesis  exaiple 

NOW,  ENTER  THE  ORDER  OF  THE  SYSTEM  (UP  TO  3). 


ENTER  THE  ELEMENTS  OP  THE  PLANT  MATRTX--A. 

„  'I  “ 

0,0 

A(2,2)  = 

-2.0 

THE  A  MATRIX  (PLANT  MATRIX! 

-i.ooE+ao  0.0 

0.0  -2.00E+00 


DO  YOU  WANT  TO  CHANGE  ANY  ELEMENT  OF  THE  MATRIX? 
n 

ENTER  THE  NUMBER  OF  CONTROL  INPUTS. 


■4  7 


ENTBR  THE  ELEHEKPS  OP  THE  DISTBIBOTION  HATSIX— B 


1.0 


0.0 


0.0 


1.0 


8(1,1)  » 
8(1,2)  » 
B(2.1)  = 
B(2,2)» 


THE  B  SATSIX  (DI3 TBIBOTIDN  MATRIX) 


I.OOE^OO  0.0 
0,0  1.00E+00 

DO  YOU  WANT  TO  CHANGE  ANY  ELEMENT  OF  THE  MATRIX? 

I 

ENTER  THE  NUMBER  OF  OBSERVABLE  OUTPUTS. 

1 

ENTER  THE  ELEMENTS  OF  THE  MEASUREMENT  MATRIX'-C. 
C(1,1)  = 


1.0 


0.0 


0.0 


C{1,2)  = 


C(2,1)  = 


2.0 


C(2,2)  = 

THE  C  MATRIX  (MSASOEEMENr  lATSIX) 


o;o'  '■  2.*3oe+oo 


i,ooE*oo  0.; 

*  1  2.  i 

DO  YOU  WANT  TO  CHANGE  ANY  ELEMENT  OF  IRE  MATRIX? 


1.0 


0.0 


C  .  0 


2.0 


YOU  HAVE  TWO  OPTIONS  AVAILABLE: 

(1)  THE  CONTROL  OPTION  FOR  SOLVING  STATE-REGULATOR 
PROBLEMS,  OR 

(2)  THE  FILTcR  OPTION  FOR  SOLVING  A  CONTINUOUS  KALMAN 
FILTER  PROBLEM. 

FOR  THE  CONTROL  OPTION,  ESTER  THE  LETTER  C. 

FOR  THE  FILTER  OPTION,  ESTER  THE  LETTER  F. 

ENTER  THE  ELEMENTS  OF  THE  CONTROL  WEIGHTING  MATRIX— R. 
R(1,  1)  = 

R(1,2)  » 


R  (2,2) 


■id 


THE  R  MATRIX  (CONTROL  HEISIPING  MATRIX) 

1,00E+00  0,0 

0.0  2. OOE+00 

DO  YOU  WANT  TO  CHANGE  ANY  ELEMENT  OF  THE  MATRIX? 
ENTER  THE  ELEMEMPS  OF  THE  SPATE  HEIGHTING  MATRIX— Q 


1.0 


1.0 


1.0 


1.0 


Q(1,1)  ’ 
Q(1,2)  * 
Q(2, 1)  = 


0(2,2) 


THE  Q  MATRIX  (STATE  NEIGHPPNG  MATRIX) 

1.00E+00  1.00E+30 

1.00S+00  1.00E+00 

DO  YOU  WANT  TO  CHANGE  ANY  ELEMENT  OF  THE  MATRIX? 


ENTER  THE  INITIAL  TIMS  FOR  THE  TRANSIENT  RESPONSE. 

0.0 

ENTER  THE  FINAL  TIME  FOR  THE  TRANSIENT  RESPONSE. 

10.0 

ENTER  THE  NUMBER  OF  POINTS  OP  THE  TRANSIENT  RESPONS 
TO  3E  PRINTED.  (<100) 

10 

M «*«*****« 

***  CONTROL  OPTION  *** 

*m****************  **** 

ENTER  THE  ELEMENTS  OF  THE  TERMINAL  BOUNDARY  VALUE 
MATRIX--?. 

0.0 

0.0 

0.0 

n  „P<2,2)  = 

0.0 

THE  P  MATRIX  (TERMINAL  BOUNDARY  VALUE  MATRIX) 


DO  YOU  WANT  TO  CHANGE  ANY  ELEMENT  OF  THE  MATRIX? 


a9 


r. 


THiNSIENI  SOLOTIDN 


TIHE  = 

1.  OOOE+01 

GAINS 

0.0 

0.  0 

0.0 

0.  0 

TIME  = 

9.  000E4-0D 

GAINS 

3,75E-01 

2.81E-01 

1.40E-01 

1.1  IE-31 

TIME  * 

8.  000E>03 

GAINS 

3. 98E-01 

2. 85E-01 

1-42E-01 

1. 12E-31 

TIME  = 

7.  000E*03 

GAINS 

4. OOS-01 

2. asE-oi 

K42S-01 

1.  12E-01 

TIME  = 

6. OOOE+03 

GAINS 

4.  OOS-01 

2. asE-oi 

1-42S-01 

1.  1 2E-3 1 

TIME  * 

5.000E4-03 

GAINS 

4. OOE-01 

2. 85S-01 

1.42E-01 

1.12E-31 

TIMS  = 

4. OOOEfOO 

GAINS 

4. OOS-01 

2. 85E-01 

1, 42S-01 

1.  12E-0 1 

TIME  = 

3. OOOEfOO 

GAINS 

4.  OOE-01 

2,  35E-01 

1.42S-01 

1,  12E-01 

TIME  = 

2.  000E«-03 

GAINS 

4.  OOE-01 

2.85E-01 

1.42E-01 

1.  12E-3  1 

TIME  = 

1. 0OOE»-O3 

GAINS 

4.  OOE-01 

2. 85E-01 

1.42S-01 

1.  12E-3  1 

TIME  = 

2. 861E-06 

GAINS 

4.  OOE-01 

2. 85E-01 

1 .42E-01 

1. 12E-31 

***********lr***  *******  ****  ******************************  *** 
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TIME  = 

2.  000E«'03 

GAINS 

2,78E-01 

2.  13E-D1 

8,53E-01 

7. 06E-01 

TIME  = 

2. 500E+0a 

GAINS 

2.79E-01 

2.  13E-t)1 

8.53E-01 

7. 06E-0t 

TIME  * 

3. OOOE+03 

GAINS 

2.79E-01 

2. 13E-D1 

8.53E-01 

7. 06E-01 

TIME  = 

3. SOOE^OO 

GAINS 

2.79E-01 

2.  13E-0  1 

8.53E-01 

7. 06E-01 

TIME  = 

4. OOOE+03 

GAINS 

2.79E-01 

2.  13E-31 

8.53E-01 

7. 06E-01 

TIME  = 

4.  5  00E+00 

GAINS 

2.79E-01 

2. 13E-01 

8.53E-01 

7,  06E-31 

TIME  * 

5.  OO0E*O0 

GAINS 

2.79E-01 

2. 1 3S-01 

8.53S-01 

7. 06E-31 

««*«**«4i****««******«4(  **4c*4i«*«***«********«***4i4i«4i*«****  ««4i 

DO  YOU  WANT  THE  TONTSOL  OPTION? 
r 

THIS  CONCLUDES  THE  RICCATI  SQUATION  PRD3RAM  (RICATI)  . 

ANALYSIS  IS  COMPLETE,  DO  YDJ  WANT  TO  RUN  LINCON  AGAIN? 
r. 

LINCON  IS  NOW  TER.IINArSD, 

The  camputer  results  far  both  optirns  presented  in 
Appendix  F  indicates  that  a  steady  state  gain  matrix  has 
beer,  achieved. 
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D.  DISCRETE  TIME  KALMAN  FILTER  (KALMAN) 


The  discrete  Kalman  filtar  gain  matrix,  G(A),  is 
calculated  by  this  program.  The  theoretical  concepts 
involved  in  the  development  of  the  computer  codes  are 
described  by  Sage  [  5].  A  briaC  development  of  the  discrete 
Kalman  filter  is  included  hara  as  an  aid  in  the  use  of 
KALMAN.  Additional  raferr'nce  materials  can  be  found  in  the 
bibliograph  y. 

Kalman  filtering  is  a  method  of  obtaining 
minimum-variance  estimates  of  signals  from  noisy 
measurements.  The  iiscrete  Kalinar.  filter  provides  state 
estimates  for  the  following  system 

i(h)  =  Is  (k-1) +4u  (k-1>  ♦PiiCh-l)  (4.D-1) 

with  the  discrete  linsar  observations 

s(k)  =  (4.D-2) 

where 

j  =  the  nx  1  state  vector  at  the  time  t(lc) 

§  =  the  nxn  nonsingular  state  transition  matrix 
r  =  the  nxr  disturbance  transition  or  distribution  matrix 


4  *  'the  nxp  control  iistributioa  matrix 

w  =  the  rxl  iisturbaace  of  sjrstaa  random  input  7ector 

5  “  the  mx1  measuramant  vector 


H  =  the  mx  1  aeasursasnt  or  observation  vector 


▼  =  the  mx 1  measurement  noise  vector 

a  =  the  px1  control  or  test  signal  vector 

k  =  the  discrete-time  index  (k  =  0,1,...) 

The  optimal  filtered  estimate  of  x(k),  denoted 
5(k|k),  is  given  by  the  recursive  relations 

l(k|k-1)  *  |5(k-l|k-1|  ♦  ^a(k-1t  (4.D-3) 


and 

i(k|k)  =  i(k|k-1)  ♦  i(k)[z  (k)-ai(k|  k-1)  ]  {4.  o-u) 

for  k  =  0,1,...,  where  x(0,0)  =  0.  The  Kalman  gain  matrix. 


G(k)  ,  is  an  nxm  matrix  which  is  specified  ae 
G(k)=P(kik-1)H^[H£(k|x-1)a^  ♦B(ki  ]-i 
£(k1  k)  =  Ci-G(k)H]P(k|  k-1) 
g(k|k-1)  =  |p(k-l|k-l)|  ♦2(k,k-1| 


(4. d-5) 
(4.  D-6) 
(4.  D-7) 
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where 


I  *  the  nxn  identity  latrix 

£<0|0>  =  E(0)^  the  initial  condition  matrix 

A  I 

j(k|k-1>  ^  the  singla-stage  prediction  of  5()c) 

j(k|)c)  *  the  filtered  estimate  of  x  (k) 

Q  *  Cs  (S  0^)  ,  t  he  ,nxn  .  co/ariance  matrix  of  the 

random  lapjt 

g  =  the  mxn  covarianca  matrix  of  the  measurament  noise 

E {w  (k)  e (k) ,  the  jeaa-square  magnitude  of  the  perturbation 
“  acceleration  matrix 

1,  InEut  Bequiram ents 

This  program  coaputas  the  recurrence  Sqs. 
(4.0-5) ,  (U.  D-5)  and  (U.0-7)  for  a  specified  number  of 
iterations  NP  and  prints  3(k.)  as  a  funcrion  of  k. 

The  required  inputs  are 

1)  the  transition  matrix  (f) 

2)  the  distribution  matrix  (P) 

31  h:{s  (k)w(ki'^  } 

4)  the  observation  matrix  (H) 

5)  the  number  of  points  •:  o  be  performed  (NP) 

and 

the  initial  condition  matrix  [P(3(0)  ]. 

5  5 


6) 


2.  g.U£Sli 


This  sectio.1  contains  a  sassior.  for  a  specific 
problem.  Commands  altered  by  the  user  are  in  lower  case. 
The  following  systam  is  co  be  tasted: 


x(k) 


.0 


0.0 


0.5  0.125 

X  (k-1)  w  (k- 1) 

1.0  *  0.5 


z(k1  =  [  1 .0  0.0]  j(k)  *■  j(k) 

R  =  [5.01 

E{w{k)w(kr}  =  ‘*•0 


£(0) 


'’1000.0 

0.  0 


0.0 


1 


1000.0 

J 


The  number  of  time  points  -^o  ba  computed  is  chosen  zc  be  20. 


iir.con 

EXECUTION  BEGINS... 

LINCON  CONSISTS  OF  THE  FOLLOWING  SUBPROGRAMS; 

BASIC  MATRIX  MANIPULAriON  -  <BASMAT> 

RATIONAL  TIME  RESPONSE  -  <RrRESP> 

STATE  7ARIABLE  FEEDBACK  -  <3rVAE> 

CONTSOLIABILITY  AND  OBSERVABILITY  -  <OBSCON> 
lUSNBESGER  OBSERVER  -  <HIEN> 

OPTIMAL  CONTROL/KALMAN  FILTERS  -  <RICAri> 

DISCRETE  TIMS  KALMAN  FILTER  -  <KALMAN> 

OPTIMAL  CONTROL  -  <OPrCON> 

PARTIAL  FRACTION  EXPANSION  -  <?RFEXP> 

ROOTS  OF  A  POLYNOMIAL  -  <RO0rS> 

TO  USE  ONE  OF  THE  SJBPROGaAMS  ENTER  THE  SAME  BETWEEN  THE 
SYMBOLS  <  >. 
kalman 

KALMAN  DETERMINES  THE  DISCRETE  KALMAN  FILTER  GAIN  MATRIX 
—  GfK)  .  FIRST  ENTER  THE  PROBLEM  IDENTIFICATION 
(*NOT  TO  EXCEED  23  CHARACTERS*), 
thesis  example 

NOW,  ENTER  THE  ORDER  DF  THE  SYSTEM  (UP  TO  8). 
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ENTER  THE  ELEMENrS  OF  THE  T8ANSITON  MArRIX— PHI 
^  PHI(1,1)« 

PHI(1,2)» 

• 

PHI(2,2)* 

1  • 

THE  PHI  MATRIX  (IHANSITION  MATRIX) 

5:88i;88 

DO  TOO  WANT  TO  CHANGE  ANY  ELEMENT  OF  THE  MATRIX? 
n 

ENTER  THE  DIMENSI3H  OF  THE  RANDOM  INPUT  VECTOR. 

1 

ENTER  THE  ELEMENTS  OF  THE  DISTRIBUTION  MATRIX— GAMMA . 
GAMHA(1  ,1)  - 
.  125 

GAMMA  (2 ,1)  * 

•  ^ 

THE  GAMMA  MATRIX  (DISTRIBUTION  MATRIX) 

DO  YOU  WANT  TO  CHANGE  ANY  ELEMENT  OF  THE  MATRIX? 
n 

ENTER  THE  ELEMENTS  OF  THE  MEA!J-52UABE  MAGNITUDE  OP  THE 
PERTURBATION  ACCSLERATION  MATRIX--W. 

THE  W  MATRIX  (MSAN-SOUARS  MAGNITUDE  DF  THE 
PERTURBATION  ACES  LERATION  MATRIX) 

U.  OOS+00 

DO  YOU  WANT  TO  CHANGE  ANY  ELEMENT  OF  THE  MATRIX? 

ENTER  THE  NUMBER  DF  OUTPUTS. 

1 

ENTER  THE  ELEMENTS  OF  THE  OBSERVATION  MATR1X--H. 

.  H(1,1)* 

1  • 

„  H(1,2)=» 

0. 

THE  H  MATRIX  (OBSERVATION  MATRIX) 

1.00E+00  0.0 
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DO  rOO  H4NT  TO  CHiNGE  ANY  ELEMENT  OP  THE  MATRIX? 

ENTER  THE  ELEMENIS  OF  THE  MEASOREMENT  NOISE  COVARIAMC£ 
MATRIX- -R. 

R(1,H» 

THE  H  MATRIX  (MEASOREMENT  NOISE  COVARIANCE  MATRIX) 
5.00E+00 

DO  TOO  WANT  TO  CHANGE  ANT  ELEMENT  OF  THE  MATRIX? 


20 


ENTER  THE  NOMBER  OF  THE  POINTS  TO  BE  PERFORMED. 

(<10  0) 

ENTER  THE  ELEMENTS  OF  THE  INITIAL  CONDITION  MATRIX--P. 


1000. 


0. 


P(1,2)  * 


0. 


P(2,  1)  * 


?(2,2)  » 

1000. 


THE  P  MATRIX  (INITIAL  CONDITION  MATRIX) 


1.  JO E*03 


03 


DO  TOO  WANT  TO  CHANGE  ANY  ELEMENT  OF  THE  MATRIX? 


K  s  0 
GAINS 
9,  95E-0  1 


0.  0 


K  =  1 

GAINS 

9.31E-01  1.92E+00 


K  -  2 

GAINS 

3.29E-Q1  9.93E-01 


K  a  3 
GAINS 

7.03E-01  6.2'4B-01 


K  a  4 
GAINS 

6.13S-01  4.52E-01 
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5 


Gkins 

5.54E-0  1 


K  «  6 

.  GMNS 
5.  18E-0  1 


K  »  7 

GAINS 
4.99E-01 


K  »  8 

GAINS 
4.  90E-01 


K  * 
GAINS 
^.87E- 


0  1 
10 


K  » 

GAI  NS 
4.  86E-01 


K  s  11 
GAINS 
4.  86E-0  1 


K  »  12 

GAINS 
4,  86B-0  1 


K  *  13 

GAINS 
4,  86E-0  1 


3.  72E-01 

3.35E-01 

3.  23E-01 

3.  23E-01 

3.  23E-01 

3.  21E-01 

3- 21E-01 

3.  213-01 

3. 213-01 


K  »  14 

GAINS 

4.  a6E-0  1 


K  a  15 
GAINS 
4. 36E-0  1 


X  a  16 
GAINS 
4.  36E-0  1 


K  a  17 
GAINS 
4.86E-  0  1 


K  a  18 
GAINS 
4.  363-  01 


3.  213-01 

3. 213-01 

3.  21S-01 

3. 21E-01 

3.  213-01 


K  *  19 

G&INS 

4.86E-01  3.21E-01 


K  »  20 

G&INS 

4.86E-01  3.21E-01 

THIS  CONCLDDES  THE  DISCRETE  lALMAN  FILTER  PROGRAM 
(KALMAN). 

ANALYSIS  IS  COMPLETE.  DO  Y03  WANT  TO  RON  LINCON  AGAIN? 
LINCON  IS  NON  TEaMINATED. 


The  computer  results  ire  prasentsd  in 


Appendix  '■ 
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E,  STME  VARIABLE  FEEDBACK  PHD3HAM  (STVARI 
Giv^n  the  linear  time- in  variant  systen 

i  (t)  »  iS(t)  ♦  fea{t)  (‘♦.E-1) 

a(t)  =  K:r(t)  -  ]  (U.E-2) 

y(tl  =  s'^X(t)  (4.S.3) 

the  following  can  be  performei  by  STVAR 

(1)  calculation  of  the  plan;  *rranfer  function,  y(s)/U(s) 

(2)  by  defining  a  ficticious  c  vector  the  internal 
transfer  function  can  be  calculated,  x;  (s) /U  (s)  ;  for 
example,  if  Xj  (s) /U  (3)  is  desired  the  c  matrix  is 
selected  with  Cj  *  1  and  all  other  c  elements  equal 
to  zero;  cr  if  X,(s)/Sj(s)  is  desired,  calculate 
X,(s)/D(s)  and  X3(s)/U(3)  and  divide  the  two 

(3)  calculation  of  the  closed-loop  transfer  function, 
Y{s)/R  (s) 

(4)  calculation  of  the  feedback  transfer  function, 
H,,  is) 

(5)  for  a  desired  closed-loop  transfer  function,  the 

controller  gain  and  feedback  coefficients  in 
addition  to  (s)  car.  be  calculated;  the 
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feedforward  gain  is  selected  so  tbat  a  zero  steady 
state  error  results  froi  a  step  input;  the  designer 
who  wishes  other  conditions  must  rescale  the  gain 
and  feedback  coefficients  appropriately;  for 
exaaple,  if  it  is  desired  to  have  K  *  1.0  but  it  is 
calculated  as  K  *  2.0  with  feedback  coefficients  of 
k  *  0.5,  k  =0.0  and  k  =1.5,  the  procedure  to 
Bodify  the  results  would  be 

2.0 

k  a - [0.5  0.0  1.5] 

“1.0 

=  [  1.0  0.0  3.0] 

All  of  the  inforaation  necessary  for  the  user  to  solve 
state  variable  feedback  problems  is  presented  in  this 
section.  However,  tie  theoretical  concepts  involved  in  the 
development  of  the  computer  codes  are  fully  described  by 
He  r  s  a  [  ^  ]  • 

The  basic  input  contains  the  problem  identification, 
matrices  A  and  b  and  the  order  of  the  plant,  n.  These  four 
inputs  are  required  regardless  of  what  open-  or  closed-loop 
calculations  are  to  be  made. 

At  this  point  STVAH  verifies  the  controllability  of  the 
system.  Three  controllability  conditions  are  possible 

(1)  complete  controllability 
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(2)  Quaerically  an  controllable 

and 

(3)  uncontrollability 

Contol lability  arises  when  ^he  controllability  matrix 

S  *  [  k  42b  ...  i'*-‘b]  (4.S-4) 

is  nonsingular,  i.e.  ,  ietg  *  D.  Even  if  the  matrix  is 
nonsingular  problems  may  arise  if  it  is  difficult  to  invert. 
To  check  this  STVAR  multiplies  the  matrix  by  its  calculated 
inverse.  The  result  should  be  the  identity  matrix.  The 
actual  matrix  product  is  compared  with  the  identity  aatrix 
to  provide  a  measure  of  uncontrollability.  If  the  maximum 
value  of  deviation  is  not  negligible,  the  plant  is 
identified  as  numerioilly  uncontrollable,  \  deviation  larger 
than  10“3  to  1 0~*  has  been  iouai  to  indicate  difficulty  by 
Melsa  [  1]. 

Osec  I-  system  is  identified  as  being 

uncontrollable,  all  open-  and  closed-loop  calculations  are 
still  performed! 

There  are  three  possible  closed-loop  options  available, 
one  for  analysis  purposes  only  and  the  other  two  for  design 
purposes.  After  choosing  the  analysis  option,  and  supplying 
STVAR  with  the  feedforward  gain  K  and  feedback  coefficient 
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matrix  Ij,  the  pcograa  calculates  the  closed- loop 
characteristic  polyaomial  aad  the  aumerator  of  the 
equivalent  feedback  transfer  function. 

The  two  design  options  are  used  to  calculate  the 
controller  gain  and  the  feedbaok  coefficients  necessary  for 
a  desired  closed-loop  characteristic  polynomial.  The 
polynomial  may  be  entered  in  either  polynomial  form  or 
factored  form. 

SilSiSS 

This  section  ccntains  a  third  order  system  for 
analysis  by  STVAR  as  presented  by  Selsa  ‘  1].  The  state 
variable  representation  of  the  plant  is  given  by 

,  '■-1.0  1.0  o.o'  ‘'o.o’ 

x(t)  *  0.0  0.0  1.0  x(t»  4-  0.0  u(‘-) 

0.0  -3.0  0.0  1.0 

L  I  L  i 

I(t)  =  [  1.0  1.0  0.0]  x(ti 

For  the  open-loop  case,  it  is  desired  to  find  the 
internal  transfer  function  X  (si /O  (s)  and  the  plant  transfer 
function  Y(s)/(J(s).  A  ficticious  matrix,  necessary  to  find 
the  internal  transfer  function,  is  then 

s  =  [0. 0  0.0  1.0] 
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In  addition,  find  tha  values  of  the  feedforward  gain  and  the 
feedbacic  coefficients  required  to  give  a  closed-loop 
transfer  function  of 


2(3+2)__ 

s3*43**6s4'4 


lincon 

EXECOTIOS  BEGINS... 

LINCON  CONSISTS  OF  THE  FOLLOWING  SOBPROGSAUS : 

BASIC  MATRIX  MANIPOLAriON  -  <BASMAT> 

RATIONAL  TIME  RESPONSE  -  <RrRESP> 

STATE  VARIABLE  FEEDBACK  -  <3rVAa> 

CONTROL UBILITX  AND  OBSERVABILITY  -  <OBSCON> 

LOENBSRGER  OBSERVER  -  <LUEN> 

OPTIMAL  CONTROL/KALMAN  FILTERS  -  <RICAri> 

DISCRETE  TIME  KALMAN  FILTER  -  <KALMAN> 

OPTIMAL  CONTROL  -  <OPTCON> 

PARTIAL  FRACTION  EXPANSION  -  <PRFEXP> 

ROOTS  OF  A  POLYNOMIAL  -  <ROOTS> 

TO  USE  ONE  OF  THE  SUBPROGRAMS  ENTER  THE  NAME  BETWEEN  THE 

SYMBOLS  <  >. 

stvar 

STVAR  DETERMINES  INTERNAL  TRANSFER  FUNCTIONS, 

THE  PLANT  TRANSFER  FUNCTION,  THE  CLOSED-LOOP  TRANSFER 
FUNCTION,  AND  THE  EQUIVALENT  FEEDBACK  TRANSFER  FUNCTION. 
IN  ADDITION,  THE  CONTROLLER  GAIN  AND  THE  FEEDBACK 
COEFFICIENTS  NECESSARY  TO  ACHIEVE  A  SPECIFIED  CLOSED- 
LOO?  TRANSFER  FUNCTION  ARE  CALCULATED. 

FIRST  ENTER  THE  PROBLEM  IDENTIFICATION 

(*NOT  TO  EXCEED  23  CHARACTERS*)  . 

thesis  exa  irple 

^  NOW,  ENTER  THE  ORDER  OF  THE  SYSTEM  (OP  TO  8). 

ENTER  THE  ELEMENTS  OF  THE  PLANT  MATRIX--A 

.  A  (1,1)* 

*  I  • 

A(1,2)» 

0  = 

Q  A(2,1)» 

0  = 

.  A  (2,3)  = 


A  (3,  1) 


0 
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.tijj 


nr' 


-3. 


A  (3,2) 
A  (3,3) 


0.  0 
0.0 


0.3  1.30E+00 

■3.00E+00  3.3 


0 

0 

1, 


DO  TOO  WiNT  TO  CHANGE  ANY  ELEMENT  OF  IHE  MATRIX? 

ENTER  THE  ELEMENIS  OF  THE  CONTROL  VECTOR — B. 

B(1)  * 


B(2) 

B(3) 


THE  B  MATRIX  (CONTROL  MATRIX) 

0.0 
0.  0 

1.  OOE  +  00 

DO  YOO  WANT  TO  CHANGE  ANY  ELEMENT  OF  THE  MATRIX? 

n 

if^*ik*Tli*****  ***********  ****  ******************************  *** 

OPEN-LOOP  CALCOLATIONS 

DENOMINATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S 
3,  OOE+00  3.00E+00  1.33E+00  1.03E+00 


THE  ROOTS  ARE 


HEAL  PART  IMA^^^PA^ 

olb  ll73E*00 

-I.OOEfOO  0.0 


0 

0 

1. 


DO  YOO  HAVE  A  FICTICIOOS  OOTPOT  VECTOR  TO  ENTER? 

ENTER  THE  ELEMENTS  OF  THE  FICTICIOUS  OUTPUT  VECTOa--C. 
C(1)* 

C(2)  = 

C(3)* 


TH 


HE  C  M. 

0.0 


MATRIX  (FICTICIOOS  OUTPUT  VECTOR) 


3.  0 
).  Q0E*Q0 

DO  YOU  WANT  TO  CHANGE  ANY  ELEMENT  OF  THE  MATRIX? 


NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S 
0.0  1.00S  +  00  1.00E  +  00 
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n 

1, 

1. 

0 

n 


P 

U. 
6 . 
u. 


THE  ROOTS  ARE  REAL  PART  IMAG.  PART 

-I.OOEfOa  0.0 

0.0  0.0 

DO  TOO  HA7S  ANOTHER  FICTICIOUS  MATRIX  TO  ENTER? 


ENTER  THE  ELEMENTS  OF  THE  TRUE  OUTPUT  7ECrOE— C. 
C(1)» 


C(2)  = 


C(3)  = 


THE  C  (OUTPUT  VECTOR) 

l!  OOE  +  OO 

0.0 

DO  YOU  WANT  TO  CHANGE  ANY  ELEMENT  OF  THE  MATRIX? 


NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S 
2.Q0E  +  00  1.00E>5o 

THE  ROOTS  ARE  REAL  PART  IMAG.  PART 

-2.0^E^^0  O.O 

THIS  PROGRAM  IS  CAPABLE  OF  PERFORMING  THREE  TYPES 
OF  CLOSED-LOOP  CALCULATIONS.  ONE  TYPE  FOR  THE  ANALYSIS 
MODE.  THE  OTHER  TWO  FOR  DESIGN. 

ENTER  ONE ‘OF  THE  FOLLOWING  MODES; 

1)  A  —  FOR  THE  ANALYSIS  MODE 

2)  P  —  FOR  THE  DESIGN  MODE  WITH  THE  UNFACTORSD 

CLOSED-LOOP  CHARACTERISTIC  POLYNOMIAL 

3)  P  --  FOR  THE  DESIGN  MODE  WITH  THE  FACTORED 

CLOSED-LOOP  CHARACTERISTIC  POLYNOMIAL 


CLOSED-LOOP  CALCULATIONS 
KEY  =  ® 

ENTER  THE  DESIRED  CLOSED-LOOP  CHARACTERISTIC  POLYNOMIAL 
COEFFICIENTS  IN  ASCENDING  POWERS  OF  S. 

YOUR  HIGHEST  ORDER  COEFFICIENT  MUST  BE  A  VALUE  OF  ONE. 
POLY  (1)  = 


POLY  (2)  = 


POLY  (3)  * 


POLY (tt)  = 


THE  NUMERATOR  OP  H-EQUIVALENT  - 
IN  ASCENDING  POWERS  OF  S 

5.00S-01  1.50E+00  1.50E*00 

THE  ROOTS  ARE  REAL  PART  IMAG.  PART 

-5.00S-01  -2.89E-01 
-5.0DE-01  2.89E-01 
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THE  FEEDBACK  COEFFICIENTS 

5.00E-01  0.3  1.50B  +  00 


THE  GAIN 


2.0  00000E«>00 


THE  CLOSED-LOOP  C HABACTERISTIC  POLYNOMIAL  - 
IN  ASCENDING  POiERS  OF  S 

4,OOE  +  O0  6.00E+00  4.33E  +  00  I.OOE^OO 


THE  ROOTS  ARE 


-I.OOEfOO  -1.00E+00 
-1.00E»30  I.OOE^OO 


MAXIMUM  NORMALIZED  ERROR  =  0.0 

DO  TOO  WANT  TO  RUN  ANOTHER  NODE  IN  SI7AR? 

E 

THIS  CONCLUDES  THE  STATE  VARIABLE  FEEDBACK  PROGRAM 
(STVAR)  . 

ANALYSIS  IS  COMPLETE.  DO  YOJ  WANT  TO  RUN  LINCON  AGAIN? 
n 

LINCON  IS  NOW  rERMINATED. 

The  conputsr  results  f^r  this  problem  is  prassntsl 
in  Appendix  H.  The  cpen-loop  portion  indicates  that  the 
system  is  controllable/  since  there  is  no  indication  of 
ur.>wOr.  troll  ability,  ^ith  the  desired  internal  transfer 
function  of 


if  I 

0  (ir’ 


s 

33 ♦s2>5s>3 


and  a  plant  transfer  function  of 


s3t.  s2+3s  +  3 


The  closed-loop  portion,  using  the  polynomial  design 


mode. 


the  feedforward  oain 


shows 


th;  t 


and  feedback 


coefficient  matrix  repaired  to  yield  the  desired  closed- loop 
transfer  function,  y(3)/H(s),  ace 

K  *  2.0 


and 


k  * 


0.5 

0.0 

1.5 


1 


I.  j 
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F.  LOENBEBGEB  OBSEBFSB  PBOGB&M  (LDEN) 

When  a  particular  closai-lcop  transfer  function  is 
desired  and  sooe  of  the  states  are  inaccessible,  LUEN  can  be 
used  to  design  a  Lusnberger  Observer.  For  example,  if  g 
measurements  are  state  variables,  an  observer  of  reduced 
dimensions  can  be  designed  to  estimate  only  those  states 
which  are  not  measured.  The  state  estimates  generated  by  an 
observer  can  be  used  as  input  information  to  a 
controller  [8].  The  block  diagram  presented  in  Fig.  4-1 
represents  the  general  form  of  tne  system  when  a  compensator 
is  placed  in  the  feedback  path. 


Fig.  4.1.  Luenberger  Observer  Block  Diagram 
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The  plant  is  characterized  ay  the 

equations 

state  and  measurement 

i(t)  *  ♦  Sa(t» 

(«.  F-1J 

2(t)  »  Ci(t) 

(4.P-2) 

2'  (t)  =  c  x(t| 

(4.  F-3) 

An  observer  can  be  designed  that  generates  an  estimate  i 

which  converges  to  the  state  x  as  time  becomes  large  [  9]. 

A  linear  controller  is  designad  as 

u(t)  =  K[E(t)  -  ic'^5(t)  ] 

(4.  F-4) 

in  which  all  states  can  be  measured 

state  with  the  estimate  yields 

.  Replacing  the  true 

u(t)  »  KCr(t)  -  !fi(t)  ] 

(4.  ?-5) 

where 

2  (hi 

(4.  F-6) 

As  tiae  increasss  ic^Z(t)  -jUI  approach  !c^x{t>.  Sines  the 
reduced-or d«r  obsarvar  astiaatss  only  the  inaccssibls 


Sta*  es , 

control 

law  contains 

measured  states 

available 

a  nd 

obser vsr 

estimates  for 

the  other  states 
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For  ease  in  using  LQEtl,  the  following  are  defined: 

X(t)  =  n-elsmant  coluan  state  vector 

u(t)  »  plant  input 

z(t)  =  q-vector  of  systea  Baasuremsnts 

z’  (t)  a  output  variable  to  be  controlled 

A  a  plant  latrix  (n  t  n) 

b  a  distribution  matrix  (q  x  m) 

[assumptions:  (1)  g  S  n  and  (2)  C  has  rank  q] 
C  a  output  aatrix 

r(t)  a  forcing  function 

K  a  feedforward  or  ccatroller  gain 

k  a  feedback  coefficiant  inatrix(n  x  m) 

h  a  observer  feedback  coefficient  marrix 

a  a  ourput  feedoack  coefficient  aatrix 

and 

2(t)  a  astimatei  stare  vector  of  x(t). 
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Froa  theorj,  the  obsacver  Is  described  by 


l(t)  »  El(t)  ♦  ♦>  SjU(t)  (<*.  F-7) 

where 

P  *  observer  eigenvalue  eatrix 

and 

g,  and  gj.  a  observer  gain  matrices. 

Z’(s)/R(s)  can  be  solved  by  the  following  procedure: 

(1)  select  Z*(3)/8(s)  and  salve  for  K  and  k  by  using  the 
STVAB  prograa 

(2)  use  the  0B3C0N  program  to  calculate  the 
observability  index;  the  observer  is  designed  to 
have  a  dimension  greater  than  or  equal  to  the 
observability  index  minus  one. 

(3)  select  the  ?  eigenvalues;  these  should  not  equal 
those  of  A,  previously  calculated  by  SIVAR. 

(4)  use  LOEN  to  calculate  G,  ,  G^,  h  and  g. 

1.  Termn^l  Session  Example 

The  example  presented  here  is  taker,  from  Desjardins, 
pp.m7  [  21.  Commands  entered  by  the  user  are  in  lower 
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esse.  However,  dae  to  the  spsrial  nature  of  the  exaiple, 
consents  have  been  added  to  the  terninal  session  for 
clarity.  Given  the  fourth  order  plant 


3.  0 

1 .0 

0.0 

0.0 

•■o.o 

0.  0 

3.0 

1.3 

0.0 

x(t>  ♦ 

0.0 

3.  0 

3.0 

0.0 

1.0 

0.0 

3.  0 

-15. 

-23. 

-9.0 

j 

1.0 

1. 

with 

z*(t)  =  [  10.  20.  0.0  0.0]  x(t> 

Because  x,  and  x^  are  the  only  aeasureable  states 


I  (t> 


1  .0  0.0  3.0  0.0 

X  (t) 

3.0  1.0  3.0  0.0  ' 

L  J 


The  closed-loop  transfer  function  to  be  achieved  is  chosen 
to  be 


_  _ 1 _ 

R(s)  ~  s*+3s3+17s2+233+20 


lincon 

EXECUTION  BEGINS... 

LINCON  CONSISTS  OF  THE  FOLLOWING  SUBPROGRAMS; 
BASIC  MATRIX  MANEPOLAIION  -  <BASMAT> 

RATIONAL  TIME  RESPONSE  -  <3rRSSP> 

STATE  VARIABLE  FEEDBACK  -  <3rVAR> 

CONTROL UBILITY  AND  OBSERVABILITY  -  <OeSCON> 
LUENBERGER  OBSERVER  -  <L0EN> 

OPTIMAL  CONTROL/KALMAN  FILTERS  -  <RICAri> 
DISCRETE  TIME  KALMAN  FILTER  -  <KALMAN> 
OPTIMAL  CONTROL  -  <OPrCON> 

PARTIAL  FRACTION  EXPANSION  -  <PRFEXP> 


7h 
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DO  TOD  M&NT  TO  CHANGE  ANT  ELEMENT  OF  THE  MATRIX? 


NOME8ATOR  COEFFICIENTS  -  IS  ASCENDING  POWERS  OF  S 
2.00E  +  01  1.00B+01 

THE  ROOTS  ABE  REAL  PARI  IHAG.  PART 

-2.00E»-00  0.0 

THIS  PROGRAM  IS  CAPABLE  OF  PERFORMING  THREE  TYPES 
OF  CLOSED-LOOP  CALCULATIONS.  ONE  TYPE  FOE  THE  ANALYS 
MODE.  THE  OTHER  TWO  FOR  DESIGN. 

ENTER  ONE  OF  THE  FOLLOWING  MODES: 

1)  A  —  FOR  THE  ANALYSIS  MODE 

2)  P  --  FOR  THE  DESIGN  MODE  WITH  THE  ONFACTORED 

CLOSED-LOOP  CHARACTERISTIC  POLYNOMIAL 

3)  F  —  FOR  THE  DESIGN  MODE  WITH  THE  FACTORED 

CLOSED-LOOP  CHARACTERISTIC  POLYNOMIAL 


IS 


20 


CLOSED-LOOP  CALC3LATIONS 
KEY  =  P 

ENTER  THE  DESIRED  CLOSED-LOOP  CHARACTERISTIC  POLYNOMIAL 
COEFFICIENTS  IN  ASCENDING  POWERS  OF  S. 

YOUR  HIGHEST  ORDER  COEFFICIENT  MOST  BE  A  VALOE  OF  ONE. 
POLY  (1)  * 


POLY  (2)  = 


POLY  (3)  * 

POLY  (4)  * 

POLY  (5)  = 

THE  NDMERATOR  OF  H-E)  DIVALENT  - 
IN  ASCENDING  POWERS  OF  S 

2.  OOE  +  01  1.30E+01  -3.00E  +  00  -3.CIOE+00 

THE  ROOTS  ARE  REAL  PART  IHAG,  PART 

-1.04E«-00  0.0 

3.21S*-00  0.0 

THE  FEEDBACK  COEFFICIENTS 

2.  OOE>01  1.30E+01  -5.00E+00  -3.00E  +  00 

THE  SAIN  *  1.  OOOOOOE  +  OO 

THE  CLOSED-LOOP  CHARACTERISTIC  POLYNOMIAL  - 
IN  ASCENDING  POWERS  OF  S 

2.00E  +  01  2.30E+01  1.70E  +  01  6.00E<-00  1.00E*00 
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THE  ROOTS  ARE  R|A^^PAJJ 

-lIOOE+OO  2l00E*00 
-2.00E*’30  0.0 

-2.00E+00  0.0 

MAXINUH  NORMALIZED  ERROR  -  0.0 

DO  TOO  HAMT  TO  RON  ANOTHER  SODE  IN  STVAR? 

THIS  CONCLUDES  THE  STATE  TARIABLE  FEEDBACK  PROGRAM 
(STVAR) . 

C(2M|ffiNT:  rasults  indicate  that  the  system 

con'wollabie,  that  the  eigenvalues  of  A  are  -3. 
-5.0,  -1.0  and  0.0,  that  the  values  of"k  are  20. 
13.0,  -6.0  and  -3.0  and  that  K  is  equal  to  unit 
the  observability  index  is  next  needed  to  design  t 
observer 

ANALYSIS  IS  COMPLETE.  DO  TOO  WANT  TO  RON  LINCON  AGAIN? 


LINCON  CON 
BASIC  M 
RATIONA 
STATE  V 
CONTROL 
LOENBER 
OPTIMAL 
DISCRET 
OPTIMAL 
PARTIAL 
ROOTS  0 
TO  OSS  ONE 
SYMBOLS  < 
Ob scon 

OBSCON 

CONTROL 

FIRST  S 
(♦NOT  T 
thesis  exa 


SISTS  OF  THE  FDLLONIN3  SUBPROGRAMS; 

ATRIX  MANIPOLATION  -  <BASMAT> 

L  TIME  RESPONSE  -  <RrRESP> 

ARIABLE  FEEDBACK  -  <3rVAR> 

U3ILITY  AND  OBSERVABILITY  -  <OBSCON> 

GER  OBSERVER  -  <LOEN> 

CONTROL/KALMAN  FILTERS  -  <RICAri> 

E  TIME  XALMAN  FILTER  -  <KALMAN> 

CONTROL  -  <OPrCON> 

FRACTION  EXPANSION  -  <PRFEXP> 

P  A  POLYNOMIAL  -  <R00rS> 

OF  THE  SUBPROGRAMS  ENTER  THE  NAME  BETWEEN  THE 

>. 


DETERMINES  THE  OBSERVABILITY  INDEX  AND 
LABILITY  OF  A  SYSTEM. 

NTER  THE  PROBLEM  IDENTIFICATION 
0  EXCEED  23  CHARACTERS*), 
sple 


NOW,  ENTER  THE  ORDER  OF  THE  SYSTEM  (U?  TO  8) 

ENTER  THE  ELEMENTS  OF  THE  PLANT  MATRIX--A. 
A(1,  1)  * 


A(1,2)  = 


A(1,3) 


A  ( 1 , 4) 


A(2,1)  = 
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0 


A  (2,  2) 


1, 

0 

0 

0 

0 

1. 

0 


A  (2,  3)  » 
A(2,4)  » 
A  (3,1)  » 
A  (3,  2)  = 
A  (3,  3)  » 
A(3,4)  - 
A  (4,1)  * 


.,^A(4,2)  = 

.23A(4,3)  = 


-9. 


A  (4,4)  = 


3:8 

0.0 

o.c 


n 

1 

0 
0 
0 
1 . 


THE  A  MATRIX  (PLANT  MATRIX) 

1.pOE+00  3.0  D.p 

0.0  1.00S+-00  3.0 

0,0  3,0  I.OOE^OO 

-t,50E+01  -2.30E  +  01  -9.00E+00 

DO  TOO  WANT  TO  CHANGE  ANY  ELEMENT  OF  THE  MATRIX? 
ENTER  THE  NUMBER  3F  COLUMNS  OF  THE  3  MATRIX. 


ENTER  THE  ELEMENTS  OF  THE  B  MATRIX. 
3(1,1)  * 


B(2,  1)  » 

3(3,  1)  * 

3(4,  1)  = 

THE  3  MATRIX 


C .  0 
0.0 
0.3 

1 . OOE+00 

DO  YOU  WANT  TO  CHANGE  ANY  ELEMENT  OF  THE  MATRIX? 
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SNTER  THE  NUMBER  OF  OUTPUTS. 


ENTER  THE  ELEMENTS 
C(1,1)  » 


OF  THE  C  MATRIX. 


C(1,2)  » 
C(1,3)  * 
C(1,4)  = 
C(2,1)  * 
C(2,2)  » 
C(2.3)  = 


C(2,4)  * 


THE  C  MATRIX 

2>00  0.5 


1 . 00E*00 

0.0 


1.0  0E+00 


DO  YOU  WANT  TO  CHANGE  ANY  ELEMENT  OF  THE  MATRIX? 


OBSERVABILITY  INOl 
THE  SYSTEM  (A.BI  ! 
DC  YOU  WANT  TO  R3N 


IS  CONTROLLABLE 

4c4c4c4i  :{c]|t4i4<*4t4[4i«<t4r*4c 

OB3CON  A3AIN? 


THIS  CONCLUDES  THE 
CONTROLLABILITY  PR 

COMMENT;  rasa 

^ - fEls  par 

observer;  aij 
(por  equal  to 
eigenvalues 


OBSERVABILITY  INDEX  AND 
OGRAM  (OB3CDN) 

Ifs  indicata  ?in  obsarvability  ini 
mizs  -tha  dasigr  of  a  second 
envaluea  of  -3.5  and  -4.0  ara  c 
•taosa  3f  the  plait)  as  obs 


ANALYSIS  IS  COMPLETE.  DO  YO'J  WANT  TO  RON  LINCON  AGAIN? 


LINCON  CONSISTS  0?  TH 
BASIC  MATRIX  MAMIP 
RATIONAL  TIME  RES? 
STATE  VARIABLE  FES 
CONTROLLABILITY  AS 
LUENBERGER  OBSERVE 
OPTIMAL  CONTROL/KA 
DISCRETE  TIME  KALM 
OPTIMAL  CONTROL  - 
PARTIAL  PRACTIDN  E 
ROOTS  OF  A  POLY NOS 
0  USE  ONE  OF  THE  SOB 
YM30LS  <  >. 
luen 


E  FOLLOWING  SUBPROGRAMS: 

ULATION  -  <BASMAT> 

ONSE  -  <RrPESP> 

DBACK  -  <3rVAR> 

D  OBSERVABILITY  -  <OBSCON> 

R  -  <LUEN> 

LMAN  FILTERS  -  <RICAri> 

AN  FILTER  -  <KALMAN> 

<OPrCON> 

XPANSION  -  <PRPEXP> 
lAL  -  <R0DTS> 

PROGRAMS  ENTER  THE  NAME  BETWEEN  TH 


H  O  I1X 


LOEN  IS  OSED  TD  DESIGN  LOENBSRGER  OBSERVERS  TO 
ACHIEVE  A  GIVEN  TLOSED-LOOP  THANSPEE  FUNCTION  WHEN  SOME 
STATE  VARIABLES  ARE  INACCESSIBLE. 


A 

MATRIX 

(PLANT  MATRIX! 

0, 

,  0 

1.00E+00 

3.3 

0. 

,0 

0, 

,  0 

0.0 

1. 33E+00 

0. 

,0 

0, 

,  0 

0.0 

0.0 

1. 

,00E*00 

0, 

,  0 

-1.50E+01 

-2.30E4-01 

-9. 

.OOE+00 

DO  YOU  WANT  TO  CHANGE  ANY  ELEMENT  OF  THE  MATRIX? 

ENTER  THE  ELEMENTS  OF  THE  DISTRIBOTION  MATRIX— B. 
3(1)* 

B(2)  = 

3(3)  = 

B(4)  = 

THE  B  MATRIX  (DISTRIBUTION  MATRIX) 

0.0 
0,  0 
0.0 

1  .  OOE  +  00 

DO  YOU  WANT  TO  CHANGE  ANY  ELEMENT  OF  THE  MATRIX? 

ENTER  THE  ELEMENTS  OF  THE  OUTPUT  MATRIX--C. 

C(1,1)  = 

C(1,2)  = 

C(1,3)  = 

C  (1,4)  = 

C(2,  1)  = 

C(2,2)  = 

C(2,3)  = 

C(2,4)  = 

THE  C  MATRIX  (OUTPUT  HATRIXl 

1.C0E  +  00  O.U  0.3  0.0 

0.0  l.OOE+00  0.  0  0.0 

DO  you  WANT  TO  CHANGE  ANY  ELEMENT  OF  THE  MATRIX? 

ENTER  THE  DESIRED  FEEDBACK  COEFFICIENTS. 

FD3K  COZF?  (1)  = 


I 
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FDBK  COEPF  (2) 


FDBK  COEFF  (3)  = 
FDBK  COEFF  (U)’* 


THE  DESIRED  FEEDBACK  COEFFICIENTS 

-6.  OOE  +  00 
-3. OOE+00 

DO  YOU  WANT  TO  CHANGE  ANY  ELEMENT  OF  THE  MATRIX? 
n 

THE  OBSERVER  EIGENVALUES  (F  MATRIX)  CAN  BE  SUPPLIED 
EITHER  IN  THE  FORM  OF  A  CHARACTERISTIC  POLYNOMIAL  OR  IN 
THE  ROOTS  OF  THAT  POLYNOMIAL. 

ENTER  EITHER  A  P  FOR  POLYNOMIAL  COEFFICIENT  FORM 
OR  AN  F  For  FACTORED  ROOT  FORM. 

f 

ENTER  THE  REAL  PART  OF  THE  ROOT. 

-3.5 

ENTER  THE  MAGNITUDE  OF  THE  IMAGINARY  ROOT. 

0 


ENTER  THE  HEAL  PART  OF  THE  SOOT. 

-4. 

ENTER  THE  MAGNITUDE  OF  THE  IMAGINARY  ROOT. 

0 

OBSERVER  EIGENVALUES  R|A^^?ARJ 

-4l0OE«-0O  O'.O 

THE  OBSERVER  CHA S ACTEHISTIC  POLYNOMIAL 
COEFFICIENTS  IN  ASCENDING  POWERS  OF  S 
1.40E+01  7.50E+00  1.00E  +  00 

*«*«4t:ti4i***  *4141  4i  4t*4>4i  4i4c4i4i  #41  **41  *#41 

THE  F  MATRIX  (OBSERVER  EIGENVALUE  MATRIX) 

-7.  50E  +  00  I.OOE+OO 
-1.40E  +  01  0.0 

THE  G1  MATRIX  (OBSERVER  GAIN  MATRIX) 

9.  55S4.01  2.  ?2S-*-0  1 

0.0  0.3 

THE  G2  MATRIX  (OBSERVER  GAIN  MATRIX) 

-3.00E4-00 
-1  .  50E  +  00 

OUTPUT  FEEDBACK  COEFFICIENTS 
2.  OOE-t-OI  8.  0  OE+00 

COMPENSATOR  FEEDBACK  COEFFICIENTS 
1 .  C0E>00  0.  3 

THIS  CONCLUDES  THE  LUEN3ERGER  OBSERVER  DESIGN  PROGRAM. 
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DO  YOO  «AHT  TO  SaM  THE  PROGRAM  AGAIN? 

ANALYSIS  IS  COMPLETE.  DO  YOJ  WANT  TO  RON  LINCON  AGAIN? 


LINCON  IS  NOW  rERMINAPED. 


The  results  of  STYAR,  OBSCON  and  LOEN  are  shown  in 
Appendix  I.  Proa  these  the  the  observer  is  given  as 


''-7.5  1 .  o' ''yj  ( t) ''  ^  *’85.  5  29.  25'"'x,  (^-)  ' 


-14.0  0.0  Y,(t) 


9.0  0.0  Xj  ( t ) 


J  L  *  J  L 


u(t)  =  1.0{r(t)  -  [20  .0  3.5  ]  '  -  [1.0  0.0]  ■»  '  { 


<*'’i 


34 


G.  OPTIMAL  CONTROL  PROGRAM  (DPTCON) 

Given  the  linear,  tiaa-invariant  systaa  represented  as 

i(t)  »  4x(t)  >  aa(t)  (4.G-1) 

OPTCON  will  ainimize  the  cost  function 

J(N)=l/2i  (N)QX  (N)*1/2  ^[X(k)  Jx  ()t)  ♦  Ru2  (k)  ]  (4.  G-2) 

k*® 

where 

X  *  state  vector 

Q  *  meaisurenent  noise  covariance  matrix  (n  x  r.< 

N  =  number  cf  time  intervals  over  which  the  sum  is  made 
R  =s  random  input  (a  scalari 
4  ®  plant  matrix  (n  x  n) 

B  =  distribution  matrix  (n  x  1) 

and 

u(t)  =  ccntrol  (a  scalar). 

The  output  of  the  program  is  the  feedback  gain  matrix 
which,  when  multiplied  by  the  state  vector,  yields  a  scalar 
ccntrol.  The  following  recursive  equations  were  derived 
using  dynamic  prograaaing, starting  at  the  terminal  time  and 
working  backwards. 

£(k)  =  (k)P(k-1)^(k|  *  2  *  aA(k)  a"'*  (k)  ,  P(3)  =  3  (4.3-3) 
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=  f  ♦  44^  (Jc),  =0 

i^(k)  =  -CAP(it-i)i3/:4^£(ic-i)4  ••  R]»  4^(0) 


0 


(4.  G-4) 


(4.  G-4) 

For  siaplicity  ia  programaiag,  the  following  terns  are 
defined; 

terminal  =  VSx"*^  (S )  Qx(H) 

trajectory  =  1/2  \  x(t)  25(^> 

w-i  ^ 
fuel  =  1/2  ^Ru2(!c) 
l|:o 

1.  Ter  ginal  Session  Example 

Given  the  aysteai  and  parameters  described  below  find 
the  discrete  steady  state  gains  for  a  sampla  of  0.1. 

.  ''o.o  i.o''  ’’o.o’ 

n  .  !c(t)  ♦  ,  ^  u(t) 

1.0  0.0  1.0 

I.  j  c  j 

'’l.O  1.0*' 

^  "  1.0  1.0 

I.  J 

F  =  1 .0 

In  addition,  run  the  program  for  a  time  interval  of 

40. 

lincon 

EXZCUTTON  BEGINS... 

3  5 


LINCON  CONSSTS  OF  THE  FOLLOMIS3  SOBPHOGRiHS ; 

BASIC  MWEIX  MANIPOLAnON  -  <BASMAT> 

RATIONAL  riME  RESPONSE  -  <arHESP> 

STATE  VARIABLE  FEEDBACK  -  <3rVAR> 

CONTROL lAB I LITT  AND  OBSERVABILITY  -  <OBSCON> 
LUENBERGER  OBSERVER  -  <LUEN> 

OPTIMAL  CONTROL/KALMAN  FILTERS  -  <RICAri> 

DISCRETE  TIME  KALMAN  FILTER  -  <KALMAN> 

OPTIMAL  CONTROL  -  <OPTCON> 

PARTIAL  FRACTION  EXPANSION  -  <PRFEXP> 

ROOTS  OF  A  POLYNOMIAL  -  <aODTS> 

TO  USE  ONE  OF  THE  SUBPROGRAMS  ENTER  THE  NAME  BETWEEN  THE 

SYMBOLS  <  >. 

optcon 


OPTCON  MINIMIZES  T 

J  (N)  =MIN  (SUM  (X  ( 

THE  OUTPUT  OF  THE 
A  TRANSPOSE,  WHICH 
YIELDS  A  SCALAR  CO 
THE  FOLLOWING  HECU 
DYNAMIC  PROGRAMHIN 
WORKING  BACKWARDS: 

(1)  AT(K)=-(DEL 

(2)  PSI{K)=P 

(3)  P(K)=PSIT(K 


HE  THE  FOLLOWING  COST  FUNCTION: 

N)  T*Q*X(N)  1-UT  (N-1)  *R*a  (N-1)  )  ) 

PROGRAM  13  THE  FEEDBACK  GAIN  MATRIX, 
WHEN  MULTIPLIED  BY  THE  STATE  VECTOR 
NTEDL 

RSIVE'eOUATIONS  were  DERIVED  USING 
G,  STARTING  AT  THE  TERMINAL  TIME  AND 

T*P(K-1)  ♦?HI/(DELT*? (K-1)  ♦DEL+R 

AT  (0)  =0 

HI*DEL*AT(K)  PSI  (0) 

)*P(K-1)  *PSI  (K)  +Q+a*A(K)*A^^(K|^ 


FIRST  ENTER  THE  PROBLEM  IDENTIFICATION 

(♦NOT  TO  EXCEED  20  CHARACTERS*).  thesis  example 

ENTER  THE  NUMBER  OF  TIME  INTERVALS  (NSTAGE)  OVER  WHICH 
THE  SUM  IS  TO  BE  MADS. 

*♦♦  NSTAGE  MUST  BE  ENTERED  IN  13  FORMAT  ♦** 
♦♦♦(I.E.,  RIGHT  JUSTIFY  TO  THREE  DIGITS  *♦♦ 


2 


1. 

0 

0 

2. 


n 

1. 


ENTER  THE  ORDER  OF  THE  SYSTEM  (UP  TO  3). 
***  ENTER  IN  II  FORMAT  *♦♦ 


ENTER  THE  ELEMENTS  OF  THE  0  MATRIX. 

***  ALL  MATRICES  ARE  ENTERED  IN  F-FORMAT  ♦** 

***  I.E.,  PUT  A  DECIMAL  POINT  AFTER  YOUR  NUMBER  *** 
Q  (^  1)  = 


Q  (1,2)  = 
Q  (2,1)  = 
Q  (2,2)  = 


THE  Q  MATRIX 
1.00E+0  0  0.0 

0.0  2.0CE>30 

DO  YOU  WANT  TO  CHANGE 


ANY  ELEMENT  OF  THE  MATRIX? 


ENTER  THE  VALUE  OF  THE  SCALAR  R 
H  = 
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ENTER  THE  SAMPLE  I 
♦  ENTER  DT  IN 


NTER7AL--0r. 

P-FORMAT 


AT  THIS  POINT  YOO  MOS' 
OPTION  A;  ENTER 

(1)  R  IS  PIMI' 
OR  IF 

(2)  R  IS  ZERO 

OPTION  B:  ENTER 

(1)  R  IS  FINI' 
OR  IF 

(2)  R  IS  ZERO 

IF  YOO  WANT  TO  REA 
BOT  NOT  THE 

HOWEVER,  IF  YOO  WA 
BOT  NOT  THE 

ENTER  THE  ELEMENTS 
.  A(1,U  = 


T  CHOOSE  ONE  OF  THE  FOLLOWING  OPTIONS 
THE  NOMBER  0  IF 

TE,  COSr*rERMINAL+TRAJECTOHY+FOEL, 

,  COST=rE3MINAL+THAJECTORY+  0 

THE  NOMBER  1  IF 
TE,  COSr*rSRMINAL+  0  FOEL, 

,  COST=TERMINAL  ♦  0  0 

D  IN  THE  \  AND  B  MATRICES, 

PHI  AND  DEL  MATRICES,  ENTER  A  0 

NT  TO  ENTER  THE  PHI  AND  DEL  MATRICES, 
A  AND  B  MATRICES,  ENTER  A  1. 

OF  THE  PLANT  MATRIX--A. 


A(1,2)  * 


1. 

0 

r. 

0 

1 . 


A  (2,  1) 

A  (2,  2)  » 

THE  A  MATRIX  (PLANT  SATRIXt 
0.0  1.00E+03 

1.00E+0  0  0.0 

DO  YOO  WANT  TO  CHANGE  ANY  ELEMENT  OF  THE  MATRIX? 

ENTER  THE  ELEMENTS  OF  THE  DISTRIBOTION  MATRIX--B. 
* 

3(2,  1)  * 


THE  B  MATRIX  (DISTRIBUTION  MATRIX) 

0.0 

1. OOE+0  0 

00  YOU  WANT  TO  CHANGE  ANY  ELEMENT  OF  THE  MATRIX? 


THE  PHI  MATRIX 

1.005003E+00  1.001657E-01 

1 .00  166'’E-01  1.005003E>-00 

THE  DEL  MATRIX 
5.00U16  UE-03 
1.00  1667E-01 


MINIMIZATION  OVER  ALL  STAGES 
•J(  1)  =AT(J)*Xp) 

AT(1)  =  -2.465-:)2 

AT(2)  =  -1.98E-01 
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2:^ 

AT  (: 

2)  • 

-3.91E-01 

oj 

AT( 

3)  » 
U  * 

AT  (; 

2)  = 

-5.79E-01 

Ih 

AT  (: 

4)  = 

1)  • 

2)  = 

-7.58E-ai 

AT  ( 

5)  * 
1)  * 

WAV 

AT  C 

2)  = 

-9.29E-51 

^  J 

AT  ( 

AT  ( 

2)  = 

-  1 .09E  +  33 

AT  ( 

f)  = 
1)  = 

AT( 

21  » 

-1.24E+33 

Oi 

AT  ( 

8)  = 
1)  * 

AT  1 

2)  = 

-1.38S<-33 

AT  ( 

lf»=“ 

AT  I 

2}  * 

-1.50E+03 

AT  ( 

10)  = 
1)  * 

AT  } 

2)  * 

-1.613+33 

AT  ( 

1  1)  = 
1)  * 

AT  ( 

2)  * 

-1 .72E+03 

AT  ( 

12)  s 
1)  * 

AT  ( 

[2)  * 

-1 .81E+33 

AT  ( 

13)  = 
* 

AT  i 

^2)  = 

-  1 .89E  +  33 

'JJ 

AT  1 

14)  = 

(1)  * 

=  AT  (J)  *Xif  J) 
-1.51E+0a 

AT  1 

[2)  * 

-1 .96E  +  03 

15)  = 
(1)  * 

'  %?  +  i3^ 

1  •  0  VD 

AT 

(2j  * 

-2.02S+03 

16)  = 
{■>)  = 

•  1  •  OO  t»  ♦  J  J 

AT 

(2)  = 

-2.07S+33 

2^ 

17)  ^ 

C)  = 

■  il.'foSlSiJ' 

AT 

hi  = 

-2.  12E  +  33 

li 

18)  ^ 
(■>)  = 

AT 

2  = 

-  2. 16E  +  33 

39 


w 


0(19)  ! 
ATf1)  » 
ATbl  “ 

iiUi  • 

Oi  21)  » 

AT(1)  = 
k1(2)  = 

0 J  22)  = 

ATfl)  = 

ATh)  » 

U(  23)  ■■ 
AT(1)  » 

ATbi  * 

Ojt  24)  = 

(1)  = 

ATb)  = 

0(  25)  = 
Af  (1)  » 
AT  (2)  » 

OJ  26)  = 

iJfll  : 

AT 


-2.20E+03 

-2.23E+03 

‘2.25E+00 

-2.27E*00 

-2.29E+a3 

-2.31E+33 

-2.32E+33 

-2.34E+33 


=  AT 


l^i 


0(  28)  = 
AT(1)  = 
AT  (2)  » 

OJ  29)  ^ 

AT(1)  = 
AT  (2)  = 

AT  (2)  = 

0(  31)  = 

AT(1)  = 

ATbi  = 

OJ  32)  ^ 

ATdl  = 


AT 


B1 


Of  3  3)  = 

AT  (1)  * 
AT  (2)  = 

Of  34)  = 

AT  (1)  = 
AT  (2)  * 

Of  35)  = 

ATfl)  * 
AT  (2  = 


-2.35E+30 

-2.35E^30 

-2.36E>03 

-2.37E+33 

-2.37E+33 

-2.38E+33 

-2.38E+33 

-2.39E+D3 

-2.39E+33 
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1 


AT  i 

(2)  » 

-2.39B+33 

AT 

hi  - 

-2.39E+03 

il 

38)  * 
(11  * 

AT 

(2)  » 

-2.UOE+35 

li 

39)  = 

(1)  = 

AT 

(2)  * 

-2.40E+33 

li 

ttO)  a 
(1)  = 

AT 

(2)  * 

-2.40E*33 

THIS  COHCLUDES  THE  OPPIMAL  EOSTROL  PHDSRAM  (OPTCON)  . 

DO  TOO  WANT  TO  HDN  THE  PR03RAM  AGAIN?  Q 

ANALYSIS  IS  COMPLETE.  DO  YDI  WANT  TO  RUN  LINCON  AGAIN? 
LINCON  IS  NOW  rERMINATED. 

The  results,  shown  in  Appendix  J,  indicate  that 


steady  state  gains  ana  achiaTsd  at  about  -2.23  and  -2.4. 
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V.  COg.SLOSiONS  AND  RECOHMENDAnONS 
A.  CONCLUSIONS 

Although  LINCON  Mas  written  primarily  as  a 
teaching/le2u:ning  tool,  it  can  still  be  guite  useful  to  the 
practicing  engineer  fDr  design  and  analysis  problems.  It  was 
written  in  modular  farm  so  that  it  could  be  easily  modified 
by  the  addition  of  subroutines. 

LINCON  has  been  extensively  tested  in  an  advanced 
optimal  estimation  course.  The  interactive  aspects  proved 
highly  successful.  Hopefully  all  the  '•bugs”  have  been 
eliminated. 

As  stated  earlier,  althouga  the  original  intent  of  this 
thesis  was  to  adapt  Desjardins'  version  of  Melsa’s  LINCON  by 
maJcing  it  interactive,  LINCON  began  to  grow  as  other 
routines  were  added  and/or  sxteisively  modified. 

OPTCON,  the  optimal  control  program  using  recursive 
equations  derived  from  dynamic  ocogramming,  is  a  new  member 
to  the  LINCON  family.  Desjardins'  KALMAN,  the  discrete  rims 
Kalman  filter  progcaa,  underwent  considerable  programming 
changes  before  achieving  i  rs  present  form. 
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B.  SSCOnMEND&TIONS 


(1)  The  source  prjgran,  LISCOS  FORTHASI  A1,  must  now  be 
passed  from  user  to  user  and  than  compilei  before  it  can  be 
used.  As  can  be  seen  in  the  terminal  session  examples,  the 
program  was  introkei  b/  typing  "iincon".  Aotually,  this  is 
an  executive  program,  LINCON  EXEC  A1,  comprising  of  the 
following  statements: 

FILEDEF  09  PRItirEP.  (RECFM  FA  LRECL  133  BLOCK  133 
LOAD  LINCON  (SIART 

The  first  statement  iafines  03  as  the  printer  and  permits  it 
to  print  out  133  characters  per  line.  The  second  statement 
invokes  the  compiled  version  of  LINCON.  It  is  recommended 
that  LINCON  be  placed  on  a  utility  disk  so  that  users  may 
link  to  it  instead  of  the  current  procedure. 

(2)  At  timjs  it  can  be  exceedingly  difficult  to 
interpret  the  tabular  output  of  some  of  the  programs.  It  is 
recommended  that  a  graphics  package  be  developed  for  RIBESE, 
PICATI,  KALMAN  and  OPTCON.  The  package  should  be  ir.terac-ive 
wirh  the  output  being  first  displayed  on  the  terminal  screen 
and  then  allowing  tne  user  to  choose  the  type  of  output, 
i.e.,  7ERSATSC,  TEKTRONIX  or  print-plct. 
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(3)  In  general,  the  prograns  are  limited  to  eighth-order 
problems.  If  the  need  should  arise  to  solve  higher  order 
systems,  this  liaioation  may  oe  removed  by  extending  the 
appropriate  dimension  statements.  The  user  must  remember  to 
alrer  the  format  statement  pertaining  to  the  output,  either 
decreasing  the  significant  figures  or  adding  a  "wrap-around" 
feature  to  overcome  the  printer  limitations  of  133 
characters  per  line. 

(4)  As  it  is  written,  LUSN  can  only  solve  for  a 
reduced-order  observer-  It  is  recommended  that  the  program 
be  modified  so  that  the  user  has  the  option  of  selecting  a 
reduced-order  observer  or  an  identity  observer. 

(5)  The  ficticious  and  real  c  matrices  of  STVAR  are 
required  to  have  the  vjtensions  1  x  n.  It  is  recommended 
that  the  program  be  modified  to  accept  a  dimension  size  of  m 
X  n. 
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BASIC  MATRIX  PRDGRAM 

PROBLEM  IDENTIPICATION:  THESIS  EXAMPLE 

****************  **************  ************  *** 


THE  A  MATRIX 

1.0  00  0  0  0E+00  0.0 

0.0  -2,000303S*00 

0.0  -5.000003E-01 

THE  DETERMINANT  OP  THE  MATRIX 
-1.5C00  OOE+00 

THE  INVERSE  3F  THE  MATRIX 
1.0  00  0  0  0E+00  0.0 

0.  0  -6.  666555  E-01 

0.0  -3.  333333E-01 


0.0 

1.  OOOOOOE+OO 
1. OOOOOOE^OO 


0.0 

6. 556665E-01 
1.333333E4-00 


**********************************  *** 

THE  MATRIX  COEFFICIENTS  OP  THE  NUMERATOR  OF  THE 
PHI(S)  MATRIX 

THE  MATRIX  COEFFICIENT  OF  S**2 
1.  0  00  0  0  0E+00  0.0 

0,0  I.OOOOOOE^OO 

0.0  0.0 


0.0 

0.0 

1. OOOOOOE+OO 


THE  MATRIX  COEFFICIENT  OF  S**1 
1.  0  00  0  0  0E+00  0.0 

0.0  -2.  OOOOODE4-00 

0.0  -5.000000S-01 


0.3 

1.000000S+00 

1.000000E+00 


THE  MATRIX  COEFFICIENT  OF  S**0 
-1.5COOOOE+00  0.0 

0.  0  1.000000B+00 

0.  0  5.  000000S-01 


0.  0 

•1.  OOOOOOE^-OO 
•2.  300000E+00 


*4141  ♦4t  ♦  **********  *  ****  ******  ************************  *** 

THE  CHARACTERISTIC  POLYNOMIAL-IN  ASCENDING  POWERS  0?  S 
1.5000  00E+00  -2.  SOOOOOE^OO  0.0  1.000000B+00 

4>4<4i4c4I*«4l*«4a>**4i4t4i4t««4>4i4i«4i*4i4i«4c«4E«4t4t4i*4c4r4i4>*4>4>4t*4<4t4i4‘4t*4t4i*4i4i4t* 

THE  EIGENVALJE3  OF  THE  A  MATRIX 
REAL  PART  IMAGINARY  PART 

8.  2  287  55S-01  0.0 

-  1,  322876E+00  0.0 

I.OOOOOOE+OO  0.0 

4^*4‘4>«*4I4I4^*  4i*4<4i4<4r*4<4<4i  ft*4‘*4>«*4I4I4l4c**4I4^4<4I4<4c4c4c*4[4I4^4I4t4^4I4‘4(*4I4c*4I*4I 

THE  ELEMENTS  OF  THE  STATE  TRANSITION  MATRIX 

THE  MATRIX  COEFFICENT  OF  EXP (  8.  228755E-0 1)  T 

0.  0  0.  0  0.0 

0.0  -6,69u317E-02  3.779633E-01 

0.0  -1. 889311 E-01  1.056‘45S+00 

THE  MATRIX  COEFFICENT  OF  EXP(  -  1 . 322875 E +  0 0)  T 
0.  0  0.0  0,0 

0.  0  1.066944  E+  00  -3.  7796a  IE-91 

0.0  1.889320E-01  -6. 69a669E-02 

THE  MATRIX  COEFFICENT  OF  EX?{  1 .  00 00 00 E <-0 0)  T 
I.OOOOOOE+OO  0.0  0.0 

0.0  -2. 264977S-36  2.351023E-06 

0.0  -5. 364413E-07  1.072884E-06 
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ON  IBM/3033{U)  NAVAL  POSTGRADUATE  SCHOOL  MONTEREY  CA 
R  M  THOMPSON  DEC  82 


UNCLASSIFIED 


F/G  9/2 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 

national  bureau  of  SIANDAROS  1963  a 


APPEgBIS  B 


P&BTIAL  PS&CIIBH  EXPANSI3!( 

PROBLEH  I DENI IP I CATION:  THESIS  EXAMPLE 

**«*«*« *41*  ************************************************ 

THE  INPOT  FONCPION  SAIN 
7.OCO0OOE+0O 

**********«******************************************>..  *** 

NOHERATOE  COEPPICIENTS  -  IN  ASCENDIN3  POWERS  OF 
2.0000E+00  2.D300E*00  1.0000E+00  1.0000E+Ow 


NOHERATOE  ROOTS 
HEAL  PART 
0.0 
0.0 

-t.OOOOOOE+0  0 


ARE 

IMA6.  PARI 
-1.  414213E4-00 
1.414213E«-00 
0.0 


*********************  ************************************* 

THE  INPOT  FONCPION  3 AIM 
7.  OOOOOOE+00 

*******************************************************  *** 

DENOMINATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S 
2.000E+00  1.000E+00  0.0  3.000S^OO  1.000E+00 


DENOMINATOR  RO 
HEAL  PART 
-3.0  18a59E+00 
-9.  1  276  04E-01 
4.6  580  98E-31 
4  ,  6  58  0  9  8E-5i 


ROOTS  ARE 

IHAG.  PART 
00  0.0 

01  0.  0 

01  -8. 308427B-01 

01  8.  308427E-01 


MOLTIPLICITY 

1 

1 

1 

1 


*********************  **********************************  *** 

RESIDUE  MATRIX  -  REAL  PARI 
5.  8  10872E+00 
3.  1  708  3  7E-01 


3.  1  708  3  7E-01 

4.  36023  IE-01 
4.  360231E-01 


0.0 

2.  061406E+00 
-2.  061406E+00 

*********************  ********************************** *** 
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USMSIS  c 


SOOTS  OF  A  POLTHOMIAL 

PaOBL^  IDENTIFICATION:  THESIS  EXAMPLE 

it*^^t*^*0m********************^i************m****** 

POLINOMIAL  COEFFICIENTS  -  IN  ASCENDIH3  POWERS  OF  S 
8.0000E-01  5.1000E+00  3.0000E  +  00  1.0000B+00 

THE  SOOTS  ABE  SEAL  PART  IMAGINARY  PAST 

-1.  1I13221E+00  1.61  6248E^5o 

-1.413221E+00  -1.616248E+30 

-1.735564E-01  0.0 

t:***  *****************¥******  ******  *** 
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BATIOM&L  TIME  SB  SPOUSE 

PROBLEM  IDEHIIPICATION:  THESIS  EXAMPLE 

**«4i****4i«  4*4E«*4t**4t**  *^m**********************1i********  *** 

THE  A  MATRIX  fPLANT  MATRIX) 

I.OCOOOOE+OO  0.0  0.0 

Q.O  -2. 000300B+00  1.000000E+00 

0.0  -5.000000E-01  1.000000B+00 


0.0  -5.000000E-01  1.000000B+00 

THE  B  MATRIX  (CONTROL  VECTOR) 

0.  0 
0.0 

1.0  00  0  0  0E+00 

THE  C  MATRIX  (OOTPUT  VECTOR) 

1.  0  000  00E4-00 
4.  OOOOOOE-01 

i.oa)oooE+oo 

THE  FEEDBACK  COEFFICIENT  VECTOR 
1.000000E+00 
1.0  000  00E+00 
0.0 

THE  CCNTROLLEH  OAIN 
3  .  2  00  0  0  0E+00 

mm**  **f***¥***********  *************  *** 

INITIAL  CONDITIONS  VECTOR-X(O) 

0.  0 
0.  0 
0.0 


'i.J35S15ol25^' 


ION  JAIN 


NUMERATOR  POLIHOMIAL  OF  R  ( S)  -  ASCEN0IN3  POWERS  OF  S 
1.  OOOOOOE4-0  0  1.  000000S<-00 

NUMERATOR  ROOTS  ARE 

REAL  PART  IMAG.  PART 

-1.  0  000  OOE+00  0.0 

DENOMINATOR  POLINOMIAL  OF  3  (S)  -  ASCENDING  POWERS  OF  S 
-2.  OOOOOOE*5o  0.  0  1.000000E+00 


DENOMINATOR  ROOTS  ARE 

REAL  PART  IMAG.  PART 

1.4  1421  3S+00  0.0 

-1.41U213E+00  0.0 

************  **********************************  *** 

THE  TIMS  RESPONSE  OF  THE  STATE  X  (T) 

VECTOR  COEFFICIENT  OF  EXP ( -5 . 00 E-0 1)  r*COS ( 1  .20E+0D ) T 
■  '  '  - 3B-01  - - 


-3. 69546 


-9.71769 


(1  .20: 

7S-01 


VECTOR  COEFFICIENT  OF  EX?1 -5  .  OOE-01 )  J’*SIN  J 1 . 20E+00 )  T 
0.0  -3. 466743e-01  -7. aO 2007E-02 


THE  VECTOR  COEFFICIENT  OP  EXP  ( 

0.0  0.0 

THE  VECTOR  COEFFICIENT  OF  EXPf 
0.0  2.6l0373E-01 

THE  VECTOR  COEFFICIENT  OF  EXPf 
0.  0  1.025083B-0T 


1. OOOOOOE+OO) T 
2.  582209E-06 

1. 4  1421  3S+00)  T 
9. 1172255-01 

-  1.414213B  +  00) T 
6. 0048543-02 


***************  ********  **********************  *** 
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OBSERVABILITY,  CONTROLLABILITY 
PROBLEM  IDENTIPI CATION:  THESIS  EKAMPLE 

iti**m*  ****************  ************************************** 


THE  A  MATRIX  (PLANT  MATRIKI 
0.0  1.  000000B«-03 

-1.000000B+00  -5.  OOOOOOE-01 

0.0  D.O 


0.0 

i.oooooaE4>oo 

I.OOOOOOE^OO 


THE  B  MATRIX 
J.gOOOOOB+00 

olo 


1. 0003 OOE^OD 
1,  OOOOOOE^OD 
3.  0 


^  JHE  C  MATRIX 

ll00O0O0E*O0 
-2.00000  0E+00 


1.  OOOOOOE>OD 
1. OOOOOOE+OD 
1. 000000E4-03 


0.0 

0.0 

0.0 


*******  ********^t*****  **********************  ************  *** 


OBSERVABILITY  INDEX  =  2 

THE  SYSTEM  {A,B|  IS  ONCONT ROLLABLE 

*******************************  *******************^f****  *** 
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OPTIMAL  CONTSOL/COMriNOOas  KALMAN  FILTER  PROGRAM 
PROBLEM  IDENTIFICATION:  THESIS  EXAMPLE 

THE  A  MATRIX  (PLANT  MATRIX! 

-1.000000E+00  3.0 

0.0  -2.  OOODOOE+00 

1.J5?oBoi."SP  (JIgTRIBDTION  MATRIX) 

0.0  1. 000300E<-03 

THE  C  MATRIX  (MS  ASOR  EMENT  MATRIX) 

1.000000E*00  3.0 

0,0  2.  000000EV03 

*«*4i*4>**«>K  m*********  *#* 

THE  a  MATRIX  (CONTROL  HEI3HTING  MATRIX) 

1,OOOOOOE>00  3.0 

0.0  2.  0003  00E«-D3 

.  MATRIX! 

1.000000E+00  1. 000300E4-03 

1.00000CE+00  1.000300E+03 

CONTROL  OPTION  ♦♦♦ 

itim********in^***** 


THE  P  MATRIX  (TERMINAL  BOUNDARY  VALUE  MATRIX) 

0.  0  3.0 

0.0  3.0 

«««  «*4i4[«4i  «  4b|i4i  41**41  *41*:^  4i4i  41*  *******  *  **********************  *** 

TRANSIENT  SOLOTIOM 
TIME  *  1.000E+31 

GAINS 


0.0  0.0 

0.0  0.0 


TIME  *  9.000E+33 

GAINS 

3.745150E-01  2.  808374E-01 

1.404037E-01  1.  n2729E-01 

TIME  =  a.OOOE+33 
GAINS 

3.982797E-01  2.8a9711E*31 

1.424856E-01  1.  117260E-01 

TIME  7.000E+33 
GAINS 

3.99705dE-01  2.  847993E-01 

1,423996E-01  1.n7371E-01 

TIME  *  6.000E+33 

GAT.NS 

3.997g5  9E-01  2.  8a7759E-01 

1. 423879E-01  K  117399E-01 

TIMS  *  5.000E+33 

GAINS 

3.  997999B-01  2.  847756fi-ai 

1.423877E-01  1.  117399S-01 


TIME  a 
GAINS 

3. 99799 9E- 0  1 
1. 423877E-01 


4.000E+33 


2.  84  7756E-01 
1. 117399E-01 


TIHE  » 

GHINS 
3.  99799  9E-01 
1.tt23877E-01 


3.000B-<-00 


2.  847756E-01 
1.  117399E-01 


TIHE  » 
GAINS 


2.000E4^30 


2.  847756E-01 
1. 117399E-01 


1.000E+DD 


TIHE  * 

GAINS 
3.  997999E-01  2. 

1.423a77E-01  1, 


847756E-01 

117399E-01 


TIHE  a  2,861E-35 
GAINS 

3.  997999B-01  2.  847756E-01 

1.423877E-01  1.  117399E-01 

******************  ******************************  *** 
********************* 

***  filter  OPTIDH 

***  ****************** 

THE  P  HATHIX  (nriTIAL  BOOSDAHT  VALOE  MATRIX) 

0.0  0.0 

0.0  0.0 

******************************************************* *** 

TRANSIENT  SOLUTION 
TIME  a  0.0 
GAINS 

0.0  0.0 

0.0  0.0 


TIME  a  5,000E-01 
GAINS 

2.393162E-01 
8.018157E-01 


2. 004539E-01 
5. 849594E-01 


TIMS  =  1.000E+00 

GAINS 

2.  70831  7B-01 
8.51  176  7E-01 


2.  1279tt1E-01 
7. 053521E-01 


TIME  a  1 . 500E+30 
GAINS 

2.  76  84  8  8E-0  1 
8.532724E-01 


2.  1331  aiE-01 
7.  O5  5743E-01 


TIMS  a  2.000E+33 
GAINS 

2.78471  5E-0  1 
8. 530354E-01 


2.  132589S-01 
7.  055843E-31 


TIME  a  2.500S+03 
GAINS 

2.789688E'-01  2.  13221  1E-01 
3.528843E-01  7.  055938E-31 


TIME  a  3.000E+03 
GAINS 

2.  79127  IE-01 
8.523297E-01 


2.  132074E-01 
7. 055972E-01 


3. 500E-*-03 


TIME  a 
GAINS 
2.79178  IE-01  2. 

3.  528109E-01  7. 


132027E-01 

055984E-01 


102 


TIME  > 


4.000E4-33 


2.79T93  8E-01  2. 

8.528078E-01  7. 

TIME  ■  4.S00B+33 

G\INS 

2.791983E-01  2. 

8.528070E-01  7. 

TIME  »  5.000E+33 

GAINS 

2.791983E-01  2. 

8.528070B-01  7. 


132320E-01 

055986E-01 


132018E-01 

0559868-01 


132318E-01 

0559868-01 


i>ii)i**^*itf^-^*tp^'tii^**nt***’*****  ******************************  *** 
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DISCBETE  TINE  ECALH&N  FILTER  PROGRAM 

PROBLEM  I DERI IF I CATION:  THESIS  EXAMPLE 

***0m*m***  Ml*********  **********************************  *** 

THE  PHI  MATRIX  (TRANSITION  MATRIX) 

1.000000E+00  5.  000D00B-01 
0.0  1. OOOOOOE*OD 

THE  GAMMA  MATRIX  (OISTRIBOTION  MATRIX) 

1.250000E-01 
5.0000  00E-01 

THE  H  MATRIX  (SE AN-SQOARE  MAGNITODE  OF  THE 

PERTURBATION  SIOELERATIDN  MATRIX) 

4.000000E+00 

THE  H  MATRIX  (OBSERVATION  MATRIX) 

1.00000CE+00  3.0 

*********************  ********************************** *** 

THE  H  MATRIX  (MEASUREMENT  NOISE  COVARIANCE  MATRIX) 
5.00  0000E+00 

THE  P  MATRIX  (INITIAL  CONDITION  MATRIX) 

1.000000E+03  3.0 

0.0  1.000300Ef03 

******* ************************************************  *** 

K  =  0 


GAINS 
9.95024  8E-01 

K  *  1 

GAINS 

9.a0772  3E-01 


K  » 

GAINS 
8.  290299E-01 


2 

!-  < 

3 


K  = 

GAINS 
7.  028207E-01 

K  =  4 

GAINS 

6.  13471  9E-01 


K  = 

GAINS 
5.543699E-01 


K  = 

GAINS 
5,  18499  IE-0  1 

K  =  7 

GAINS 

4.99227  8E-01 

K  *  8 

GAINS 

4.  90446  SE-OI 


3.  0 

1.  923759E<-0D 

9.  98075  IB-01 

5.  237227E-01 

4. 522047E-01 

3.  716316E-01 

3. 3625 13E-01 

3. 232393E-01 

3. 201340E-31 


K  *  9 

GAINS 

4.  872549E-01 


3.  202326E-01 
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K  *  10 

GAIHS 

4.  8644  00E-01 

K  =  11 

GAINS 

4.863475E-01 

K  *  12 

GAINS 

4.  86345  6E- 01 

K  *  13 

GAINS 

4.86294  4E-01 

K  =  14 

GAINS 

4.  86216  eE-01 

K  =  15 

GAINS 

4.  861480E-01 

K  =  16 

GAINS 

4  .  86  1  0  3  5E-01 

K  *  17 

GAINS 

4.860806E-01 

K  =  18 

GAINS 

4,  860712E-01 

K  =  19 

GAINS 

4,a60685E-01 

K  *  20 

GAINS 

4.  860680E-01 


3. 207918E-01 

3.  2103  92E-01 

3.  2093  99B-01 

3.  20  8563E-01 

3.  2072 75B-01 

3.  2064 94E-01 

3. 206143S-D1 

3.  206035E-01 
3. 206024S-01 

3. 206034E-01 

3. 206041E-D1 


«4i  *4^  «4t4i4i  4i  *4i4c  *  4^4^  ***  4c  iK*  41  4>4c4t 
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SlilE  VABIABLS  FEEDBACK  P33GRAN 

PROBLEH  IDENIIFICATIOH;  THESIS  EXAMPLE 

4i*«**4>4i4t*4i  ******** 4i*4>  IHHm^Htm***^*************  *** 

THE  A  MATRIX  (PLANT  HATRIXI 
-1.0C00O0E+0O  I.OOOOODE+OO  0.3 

0.0  0.0  1.000000E+00 

0.0  -3.00000DS+00  0.0 

THE  B  MATRIX  (CONTROL  MATRIX) 


THE  J  MATRIX  (CONTROL  MATRIX) 

olo 

1.  OOOOOOE+00 

***  *tt*m**************  ***************************  *******  *** 

OPEN- LOOP  CALCOLATIONS 

DENOMINATOR  COEFFICIENTS  -  IN  ASCENDING  POSERS  OF  S 
3.0000E+00  3.0000E4-03  1.0000B  +  00  1.0000E+00 

THE  ROOTS  ARE  REAL  PART  IMAGINARY  PART 

0.3  -1.732051E«-00 

0.3  1.732051E+00 

-1.D00000E+00  0.0 

^HfiUm***  ************************  *************  ***********  ** 

THE  C  MATRIX  (FICTICIOUS  OOTPOT  VECTOR) 

0.  0 
0.0 

1.  0  000  OOE+00 

NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  POSERS  OF  S 
0.0  1.000303E+00  1.000000E+00 

THE  ROOTS  ARE  REAL  PART  IMAGINARY  PA.  ' 

-1.300000E+00  0.0 

3.3  0.0 

*  **4i4ii|i4i  *#**««*  4>  ****  m*****  ************  ************  *** 

THE  C  MATRIX  (03TPUP  VECTOR) 
t.OOOOOOE+00 

1.  OOOOOOE+00 

0.0 

NUMERATOR  COEFFICIENTS  -  IN  ASCENDING  POWERS  OF  S 

2.  000000E+05  1.000003E4-00 

THE  ROOTS  ARE  REAL  PART  IMAGINARY  PART 


THE  ROOTS  ARE  REAL  PART  IMAGINARY  PART 

-2. 300000E+0C  0.0 

***inm**  *******  ************************************* 

CLOSED-LOOP  CALCULATIONS 

KEY  »  P  ***** 

THE  NUMERATOR  OF  H-EQUIVALENT  - 
IN  ASCENDING  POWERS  OF  S 

5  .  0  00  0  0  0E-3  1  1.500000E  +  00  1.500000E  +  00 


THE  ROOTS  ARE 


REAL  PART  IMAG.  PART 

-5, 30a000E-01  -2.886756E-01 

-5, 300000E-01  2.886756E-01 


THE  F! 
5.  0< 


ACK  COEFFICIENTS 
OE-01  0.3 


1. 530000E+00 


THE  G  AIN  » 


2. 000300B+03 


THE  CLOSED-LOOP  CHARACTERISTIC  POLYNOMIAL  - 
IN  ASCENDING  POWERS  OF  3 

4.  0000E*00  6.0000E»>03  i».OOOOE*00 


1.0003E+00 
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THE  SCOTS  ARE 


MAXI  a  OH  NOHMALIZED 


REAL  PART  IHAG. 
-2.D00000E+00 
-1.000000E+00 
-1.003000E+00 

ERROR  =  0.0 


n? 


PART 

0.0 

-1.000000E+00 

I.OOOOOOE+OO 


SI&rB  VARIABLE  P 
PROBLEM  IDEMTIPI 

*  *«««««  4I  **  Ml  *4i  *«**«*  « 

"“S.8  ”"??iSoi5ii 

0.0  0.0 

0.0  0.0 

0.0  -1.5000 

THE  B  MATRIX  (00 
0.  0 
0.0 
0.  0 

1.  OOOOOOE+00 

************** 

OPEN- LOOP  CALCOL 
DENOMINATOR  COEP 
0.0  1.500S+0 

THE  ROOTS  ARE 


********************* 

THE  C  MATRIX  HU 
2.  OOOOOOE+OT 
1-  OOOOOOE+01 
0.0 
0.  0 

NUMERATOR  COSPFi: 
2.  OOOOOOE+01 


EEDBACK  P30GRAN 
CATION:  THESIS  EXAMPLE 

4>*4i4i *****#*« ****41411^** **««*4i4i«****«i  *** 

ANT  MATRIX) 

OE+00  3.0  0.0 

1. 000000E4-00  0.0 

0.0  1.000000E+00 

00E«-01  -2.300000E+01  -9.  OOOOOOE+00 

NTROL  MATRIX) 


41***  *4114**  4141**41**4141  #**41  *** 

ATIONS 

PICIENTS  -  IN  ASCENDING  POBERS  OP  S 
1  2.300Bf01  O.OOOEfOO  1.000S-*-30 


REAL  PART 
-3.300000S+00 
-4. 999999E4'00 
-9.999999E-01 
0.0 


IMAGINARY  PART 
0.0 
0.0 
0.0 
0.0 


**********  ************************  **** 

TPUr  VECTOR) 


THE  ROOTS  APE 

***  *******  *********** 

CLOSED-LOOP  CALC 
KEY  =  P  ***** 
THE  NUMERATOR  OF 
2. 0  0OE4-01  1 


REAL  PART  IMAGINARY  PART 

-2.000000E+00  0.0 

**********************************  «* 

ULATIONS 

H-EQUIVALENT-IN  ASCENDING  POWERS  OF  S 
.300S4-01  -6.000S4-00  -3,OOOS+O0 


THE  HOOTS  ARE 


3.206289E+00  0.0 


THE  FEEDBACK  COE 
2.  0003*01  1 

THE  GAIN  =  1. 

THE  CLOSED-LOOP 
IN  ASCS 
2.a00E4-0  1  2.800E4- 

THE  ROOTS  ARE 


MAXIMUM  NOEMALIS 


FFICIENTS 
.  300E4-01 

0003  00E4-0) 


-6. OOOE+30 


-3.000E+00 


CHARACTERISTIC  POLYNOMIAL 
NDING  POWERS  OF  S 
01  1.700S4-01  6.000E+00 


1 ,OOOE*00 


REAL  PART  IMAG.  PART 

-1.  OOOOOOE4-00  -2.COOOOOE4-00 
-1.  OOOOOOE+OO  2.  000000E4-00 

-2.300484S+00  0.0 

-1.  9995  16E+00  0.0 

ED  ERROR  =  0.0 


OBSEH VABILITY,  C ONTROLLA3t LI TY 

PROBLEM  IDENTIFICATION;  THESIS  EXAMPLE 
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******miL)tt****00*m**t*  **********  ***********************^1  *** 


o.r  *  “"■'?!oi5SS8iJ8’^'”!,.o 

0.0  0.0  1.  000000E4-00 

0.0  0.0  0.0 

0.0  -1.5000  00E«-01  -2.300000E+01 

Q  JHE  B  MATRIX 

olo 

0.0 

I.OOOOOOE^OO 


0.0 

0.0 

1.000000E+00 

-9.000000E+00 


1.  J?5oSce+S?^*  0.  0  0.0  0.0 

0.0  1.  000000B+00  0.0  0.0 

*********************  ********************************** *** 
OBSEHVABILITX  I5IDEX  »  3 

THE  SYSTEM  (A,B|  IS  CONTROLLABLE 
*********************  ********************************** *** 


LOENBERSER  OBSERVER  DESISH  PROGRAM 
PROBLEM  IDENTIFICATION:  THESIS  EXAMPLE 

**********  ***********  ********************************** *** 

THE  A  MATRIX  (PLANT  MATRIX) 

0.0  1.0(J03  00E*00  0.0  0.0 

0.0  0.0  0000002*00  0.0 
0.  0  0.0  0.  0  1.000000S+00 

0.0  -1.5000  00E4-01  -2.  300000E+01  -9. 000000  E*00 

THE  B  MATRIX  (01 STRIBOTION  MATRIX) 

0  •  0 
0.0 
0.0 

1.  OCOOOOE+00 


THE  C  MATRIX 
1.  OCOOOOE^O 
0.  0 


pOTPOT  MATRIX) 

0  0.0 

1.000000E  +  00 


[■HE  DESIRED  FEEDBACK  COEFFICIENTS 
2  .  0  00  0  0  0E+01 
1.  300000E+01 
-6.  OCOOOOE+00 
-3.  OCOOOOE+00 


OBSERVER  EIGENVALOES 


REAL  PART  IMAG.  PART 

-3.500000E+00  0.0 

-4. OOOOOOE+00  0.0 


OBSERVER  CHARACTERISTIC  POLYNOMIAL 
COEFilCIENTS  IN  ASCENDING  POWERS  OF  S 

1.4C0000E+01  7.500000S<-00  0000002*00 

THE  F  MATRIX  (OBSERVER  SI3SNVALDE  MATRIX) 

-7.  5  COO  OOE+00  I.OOOOOOEi-00 

-1.  4  00  0  0  0E+01 


-7.  5  COO  OOE+00  I.OOOOOOEi-OO 

-1.  4  00  0  0  0E+01  0.0 

THE  G1  MATRIX  (OBSERVER  3AIN  MATRIX) 
8.  549997E+0  1  2.  924994E*01 

0. 0  0. 0 

THE  G2  MATRIX  (OBSERVER  GAIN  MATRIX) 
-2.  999997E+00 
-1.  5000  01E+00 
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OOTPOT  FEEDBACK  COEFFICIESTS 
2  .  0  00  0  0  0E>01  8.500004S4-00 

COMPENSATOR  FEEDBACK  COEFFICIENTS 
1.  OQOOOOE4-00  0.0 


113 


i££Si>0£S  J 


Shif 

AT  (21 

ATb) 

0  jt  28 

mi\ 

ATb) 

0  (  30 
ATM) 
AT^2) 

or  31 
AT(1) 
AT|2f 

AT(1) 
AT  12) 

0  (  33 
AT(1) 
ATU) 

a(  34 
AT(1) 
AT(2{ 

tl(i 

a  (  36 
AT  (1) 
AT^) 

0  (  37 
AT(1) 
Albi 

0  (  38 
AT(1) 
At|2) 

Of  39 
ATM) 
AT  2 


»  -2.35E«-50 

=  -2.3SE+3  0 

=• 

a  -2.36E«-30 

a  -2.37E«'30 

a  -2.37E+00 

a  -2.33E«-30 

a  -2,38E«'30 

a  -2.39SfaO 


-2.39S«- 


a  -2.39E*30 

a  -2.39E«-30 

a  AT(JI*X(J) 
a  -2.22S«-30 
a  -2.40EV30 

a  -2.tt0E«-30 


liAV  - 

AT]2)  a  -2.40E«-30 

4i4>*4i4i*a*a*«4tik*4i**«*4^4i  *41  4i  *  «]|i  4t* ##4i4i  #  41*41  *  4t4t4i 
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